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The Global Modeling and Simulation Branch (GMSB) of the Laboratory for 
Atmospheric Sciences (GLAS) Is engaged In general circulation modeling studies 
related to global atmospheric and oceanographic research. The research activities 
are organized Into two disciplines: Global Weather/Observing Systems and Climate/ 
Ocean-Air Interactions. 

The Global Weather activities are grouped In four areas: 1) Analysis and 

Forecast Studies, 2) Satellite Temperature and Wind Retrievals, 3) Analysis and 
Model Development, 4) Atmospheric Dynamics and Diagnostic Studies. 

The GLAS Analysis/Forecast/Retrleval System was applied to the FGGE Special 
Observing periods 1 and 2 (5 Jan - 5 March and 1 June - 31 July 1979) respectively. 
The resulting analyses have already been used In a large number of theoretical 
studies of atmospheric dynamics, forecast Impact studies and development of new 
or Improved algorithms for the utilization of satellite data. The results 
obtained with the GLAS temperature retrieval system are particularly noteworthy: 
not only are the atmospheric temperature soundings quite accurate, but as a by- 
product, remarkably accurate determinations of land and sea surface temperatures, 
cloud cover, height and temperatures, and snow and Ice extent have been obtained 
from operational TIROS-N sounders. 

The focus of the climate research Is to understand the variability of the 
climate on the monthly to seasonal time scale. In support of this goal, studies 
are carried out In the three areas of: 1) Data Analysis, 2) Climate Modeling, 

and 3) Sensitivity Experiments. The analysis of observational data to formulate 
and validate hypotheses for short-term climate change has been and continues to 
be an integral part of the climate research effort. Extensive studies have 
been completed which establish the physical basis for the prediction of monthly 
and seasonal climate changes through the use of anomalies In sea surface 
temperature and other space observed boundary conditions. In conjunction with 
these studies, existing general circulation/climate models are being extended 
and improved substantially. A 24-month simulation carried out In the presence 
of a seasonal cycle Indicates that the current model can simulate the winter to 
summer climate change. Currently, the models are being extended to Include 
much or the stratosphere. 

Ocean/air Interaction studies are concentrated on the development of models 
for the prediction of upper ocean cunrants, temperatures, sea state, mixed- 
layer depths, and upwelling zones,* and',dn studies of the interaction of the 
atmospheric and oceanib circulation sy.stems on time scales of a month or more. 

This research review contains a synopsis of extended abstracts In the 
Global Weather/Observing System studies, and Climate/Ocean -Air Interaction 
studies conducted at the Modeling and Simulation Branch for 1982-83 (Sections 
2 and 3). Abstracts of the Summer Lecture Series, jointly sponsored by GMSB 
and the University of Maryland, are present In Section IV. Section V lists the 
recent publications by the GMSB and visiting scientists. 
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EMPIRICAL CORRECTION OF THE LARGE SCALE FORECAST BIAS OF A GLOBAL NMP MODEL 

R. N. Hoffman and J. Firestone 


The most satisfying solution to the problem of large scale forecast errors 
would be to find and correct their root causes. The least satisfying solution 
would be to remove the observed bias at each forecast time, as an a posteriori 
statistical filter as suggested i^y Leith (1974a). We would like to propose, 
as an intermediate solution, th# empirical correction of the large scale fore- 
cast bias within a NWP model. 

The general topic of the empirical correction of NWP models has been 
discussed extensively by Leith (1974b, 1978) and by Faller and his collabora- 
tors (Faller and Lee, 1975; Faller and Schemm, 1977; Schemm and Faller, 1977; 
Schemm et. al. 1981). Leith has theoretically solved the complete problem of 
preserving the climate mean and single-time covariance by adding constant and 
linear terms to the equations of motion. Leith's solution is not practical 
because estimates of the climate statistics are based on relatively small samples, 
and Leith's plan involves estimating 0 (n 2) parameters, where N is the number 
of model variables. If we simplify Leith's scheme to only eliminate the model's 
bias we need only estimate 0(N) parameters. But even N is much greater than 
our effective sample size. 

Both Leith and Faller suggest making corrections locally in grid point 
space because the corrections are presumably needed to account for subgrid 
scale processes which have been parameterized and truncation errors which 
primarily effect the small scales. Unfortunately, from the practical viewpoint 
the small scales are not well observed and without a good estimate of the 
error we cannot correct it. Therefore we plan to make corrections in the 
normal mode representation of the model, correcting only the largest spatial 
scales and longest time scales. There are several advantages to this approach: 

(a) The number of parameters to be estimated is greatly reduced, avoiding 
in part the sampling problem. 

(b) By excluding from consideration those normal modes of essentially 
computational or gravity wave nature we avoid estimating those parameters 
which we are unable to observe. 

(c) Restricting the corrections to Rossby modes is also an advantage 
when applying the corrections. If we correct the governing equations then we 
avoid adding a forcing term to the gravity wave modes. If we correct the 
model state at intervals the shocks to the system do not directly change the 
fast modes, although they may excite gravity waves indirectly. 

(d) A good deal of the bias appears to be large scale. 

An interesting connection between empirical bias corrections and the 
specification of topography is noted. Since topography appears in the hori- 
zontal momentum equation at all levels as a single constant the "barotropic" 
component of the corrections to the horizontal momentum equations may be inter- 
preted in terms of a topographic correction. Interestingly it turns out that 
this "barotropic" component of the wind errors is insignificant compared to 
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the wind errors at any one sigma level. 

We have been trying to Identify and plan to remove the model bias by 
examining the perceived forecast error observed during the assimilation cycle, 
1.(S. the analysis minus the 6 hour first guess, of the 6LAS Analysis/Forecast 
system (Halem et al ., 1982) for the FGGE SOP-1, The particular analysis cycle 
we are examining (2453) Is a "full FGGE" experiment, l.e., all data sources are 
Included (Baker et al., 1981a). This experiment makes use of the thoroughly 
edited FGGE II-b"Hata set (Baker et al .. j981b). 

The results to date are Intriguing; some error statistics are presented 
In the figures. These results are for 5 January through 5 March. Many of the 
small Isolated features In these figures are proba&ly due to station biases, 
but quite a few are associated with geographic factors which plausibly could 
cause some of the parameterizatlons used In the model to be In error. In 
Fig. 1 we see that surface pressure Is forecast high over Greenland, the 
Rockies and east central Alaska. In the wind field at sigma level 8 a 
distinctive southeastly flow Is seen In the error field over the North 
Atlantic (Fig. 2a), while at sigma level 3 a very well defined antlcyclonlc 
flow is seen In the error field off the West Coast (Fig. 2b). In Fig. 3 
we see that the model Is too warm over Greenland and too cold around the 
southern boundaries of Greenland. This feature Is seen from sigma level 9 
up to sigma level 5 with lltle change In shape. Figure 4 shows the rms 
temperature errors over all sigma levels and all times In the sample. The 
feature centered over the foothills of the Himalayas indicates the model Is 
regularly making poor 6 h forecasts here. The greatest part of this feature 
Is due to errors In the lowest layers. 

We are currently working on methods of filtering the observed bias; 
we are hopeful that we can separate the model bias from the station biases. 
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CASE STUDIES OF COASTAL CYCLOGENESIS 
F;« Atlas and R. Rosenberg 


A series of forecast experl laents was conducted to assess the accuracy of 
the GLAS Model's prediction of two cases of coastal cyclogenesis and to determine 
the importance of large-scale dynamical processes and diabatic heating to the 
cyclone developmeiit* The model used in the assimilation and to generate the 
forecasts is the fourth order global atmospheric model described by Kalnay- 
Rivas et al. (1977), Kalnay-Rivas and Holtsma (1979) and Haleni et aU (1982), 
with nine vertical layers, equally spaced in sigma, and a coarse lioHkontal 
resolution of 4® latitude by 5® longitude, which is compensated by the use of 
accurate horizontal differences. 

The first case studied was the "President's Day Cyclone" which developed 
along the east coast of the United States on February 18-19, 1979, and produced 
very heavy snowfall from Virginia to southeast New York. This storm was 
significant because of the severity of the weather it produced and the failure 
of the operational models in use at the National Meteorological Center (NMC) to 

adequately predict the intensity of cyclogenesis. For example, the Limited- 

Area Fine Mesh (LFM-II) model's 12, 24, 36 and 48 h forecasts which verified at 
1200 GMT 19 February 1979 all displayed serious errors in the prediction of 
this cyclone and associated heavy precipitation. The 12 h forecast predicted 
the low to be too weak and south of its observed position; the 24 h forecast 

predicted a imich weaker low; and the 36 h forecast indicated only a trough with 

no closed cyclone center. 

The GLAS model 36 h forecast from the GLAS analysis at 0000 GMT 18 February 
correctly predicted intense coastal cyclogenesis and heavy precipitation (Fig. la). 
When this forecast was repeated without surface heat and moisture fluxes, only 
an inverted trough developed along the east coast (Fig. lb), while an experiiiHjnt 
without surface heat and moisture fluxes failed to predict any cyclonic develop- 
ment (Fig. Ic). An extended-range forecast from 0000 GMT 16 February as well 
as a forecast from the NMC analysis as a forecast from the NMC analysis at 0000 
GMT 18 February interpolated to the GLAS grid (not shown) predicted weak coastal 
low development. 

Detailed examination of these forecasts (Atlas anu Rosenberg, 1982) showed 
that diabatic heating resulting from oceanic fluxes significantly contributed 
to the generation of low level cyclonic vorticity and the intensification and 
slow rate iwvement of an upper level ridge over the western Atlantic. As an 
upper level short-wave trough approached this ridge, diabatic heating associated 
with the release of latent heat intensified, and the gradient of vorticity, 
vorticity advection, upper level divergence, and upward vertical motion (Fig. 2) 
in advance of the trough were greatly increased, providing strong large-scale 
forcing for the surface cyclogenesis. 

The second case studied was the intense cyclone which developed along the 
east coast of the United States on 6-7 February 1978. This storm, which also 
produced record-breaking snowfall accumulations, was accurately predicted by 
NMG's operational forecast models (Brown and Olson, 1978; Cressman, 1978). 



P1g. 3a shows the 6LAS model 36 h forecast generated from the NMC analysis 
at 1200 GMT February 1978 while the corresponding forecasts with surface n;* 'sture 
fluxes withheld and surface heat and moisture fluxes withheld are shown In 
Figs. 3b and c, respectively. Comparison of these forecasts shows some weakening 
of the predicted cyclone Intensity resulting from the exclusion of surface 
fluxes. However, the effect of these fluxes Is significantly smaller than In 
the President’s Day cyclogenesis, 


REFERENCES 

Atlas, R., and R. Rosenberg, 19825 Numerical prediction of the Mid-Atlantic 
states cyclone of 18-19 February 1979. NASA Tech- Memo. 83992 , 53 pp. 

Brown, H. E., and 0. A. Olsen, 19785 Performance of NMC In forecasting a record 
breaking winter storm, 6-7 February 1978. Bull. Amer. Meteor. Soc., 59, 
562-575. 

Cressman, 6. P., 1978s A notable forecasting achievement. Bull. Amer. Meteor. 
lOCi, 59, 370. ■ 

Halem, M., E. Kalnay-Rivas, W. E. Baker, and R. Atlas, 1982: An assessment of 
the FGGE satellite observing system during SOP-X. Accepted by Bull. Amer. 
Meteor. Soc. 

Kalnay-Rivas, E,, A. Bayliss, and 0, Storch, 1977: The 4th order 6ISS model of 
the global atmosphere. Beitr. Phys. Atmos. , 50 , 299-311. 
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Flq. 1. GLAS model 36h sea level pressure (solid lines) and 1000-5 
(dashed lines) forecasts valid 1200 GMT 19 February 1979 w 
(a), no surface moisture flux (b), and no surface heat or 
(c). 
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THE EFFECT OF DATA AND DATA ANALYSIS TECHNIQUES ON NUMERICAL WEATHER FORECASTS 
W. E. Baker, R. Atlas, M. Halem, and J. Susskind 


In this study we examine the sensitivity of forecast skill to Individual 
components of the First GARP (Global Atmospheric Research Program) Global 
Experiment database as well as to some modifications In the data analysis 
techniques. Several short assimilation experiments (0000 GMT 18 January 1979 
through 0000 GMT 21 January 1979) are performed In order to test the effects of 
each database or analysis change. Forecasts are then generated from the Initial 
conditions provided by these experiments. The 0000 GMT 21 January 1979 case Is 
chosen for a detailed Investigation because of the poor forecast skill obtained 
earlier over North America for that particular case. Specifically, we conduct 
experiments to test the sensitivity of forecast skill to: 

1) the Individual observing system components, 

2) temperature data obtained with different satellite retrieval 
methods, and 

3) the method of vertical Interpolation between the mandatory pressure 
analysis levels and the model sigma levels. 

For the one case examined, NESS (National Environmental Satellite Service) 
TIROS-N Infrared land retrievals are found to degrade the forecast, while 
TIROS-N retrievals produced by GLAS with a physical Inversion method, as part 
of an analysis forecast cycle, result In an Improved forecast. The use of 
oceanic VTPR (Vertical Temperature Profile Radiometer) satellite retrievals 
also results In an Improved forecast over North America. The forecast Is found 
to be sensitive to the method of vertical Interpolation between the mandatory 
pressure analysis levels and the model sigma levels. 

These results are Illustrated for the $i score In Figs. 1 and 2 for 
the sea level pressure and 500 mb height forecasts verified against the ECMWF 
(European Centre for Medium Range Weather Forecasts) analysis over North America. 
In the control assimilation (Experiment 3), (rawinsondes, pibals, aircraft, 
surface and ship data are utilized). In Experiments 4 and 5 NESS TIROS-N 
retrievals are added to the control globally and only over oceans, respectively. 
VTPR data are similarly utilized In Experiment 7 over the oceans. GLAS retrievals 
are added to the control In Experiments 9 and 10 In a non-InteractIve mode and 
In an Interactive mode, respectively. In the former, the satellite temperature 
retrievals utilize a first guess which Is Independent of the analysis/forecast 
cycle In which the retrievals are used. In the latter, satellite temperature 
soundings are generated during the course of the analysis/forecast cycle. The 
model 6 h forecast provides the first guess for the retrievals as In the non- 
Interactlve cycle, but In the Interactive cycle, the current model 6 h forecast 
provides the first guess for the retrievals rather than the 6 h forecast from 
the control experiment performed earlier. 

In the experiments just described, differences between the model first 
guess and the analyzed fields are Interpolated from pressure to sigma rather 
than the analyzed fields themselves In order to reduce Interpolation errors 
during the assimilation cycle. In order to test this change, the analyzed 
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wind and temperature fields are interpolated in Experiment 11. Otherwise, 
Experiment 11 is identical to the control. As may be seen In Fig. 2, Experiment 
11 is poorer than Experiment 3. 














HdDBL SENSITIVITY TO LOW-LEVEL WIND SPECIFICATION 
R. Atlas and A. Pursch 




For Seasat-A scatterometer (SASS) data to have a significant impact on 
numerical weather forecasting, substantial differences between the atmospheric 
state analyses produced with and without SASS data must occur, and the 
assimilation/forecast model used must be sensitive to the low 1evol wind 
modifications. Assimilation of scatterometer data can work In two ways to 
modify the analysis. First, it can modify the wind field directly and the mass 
field indirectly through geostrophic adjustments. Second, these differences 
may be amplified at the next analysis time due to their effects on the first 
guess fields used for that analysis. Thus, the differences in objective analyses 
produced by the assimilation of scatterometer data can grow until a statistical 
equilibrium value is reached. This depends on the model's climatology and on 
the availability of conventional data. 

Preliminary assimilation/forecast experiments with SASS data utilizing the 
coarse 4® latitude x 5° longitude GLAS Fourth-Order Forecast Model (Kalnay- 
Rivas et al ., 1977 and Kalnay-Rivas and Hoitsma (1979) have shown a generally 
small effect of the SASS data on both analyses and forecasts. This may be due 
to the quality of the dealiased SASS winds, the assimilation methodology used 
in the preliminary Impact experiments, model iftsensitivity to the SASS winds 
resulting from coarse resolution and the lack of a detailed planetary boundary 
layer formulation, or the SASS winds may be largely redundant to existing 
conventional observations and satellite temperature soundings. 

Experiments have been conducted to investigate the sensitivity of the GLAS 
Fourth Order Forecast Model to low level wind specification. In the experiments, 
model forecasts were generated from initial conditions in which the correct 
1000 mb or 1000 and 850 mb wind fields were replaced by the corresponding fields 
from 24 hours earlier. This has the effect of inserting large spatially coherent 
errors into the initial low level wind fields without disturbing the low level 
mass field or the upper level mass and wind fields, and allows the model 
sensitivity to these errors to be examined. 

Fig. 1 shows the initial sea level pressure field for a forecast from 0000 
GMT 9 September 1978, the 1000 mb wind field valid at this time, and the 1000 
mb wind field valid at 0000 GMT 8 September. Fig. 2 shows the 36 hour forecasts 
from 0000 GMT 9 September which utilized these 1000 mb wind fields as well as a 
third forecast in which both the 1000 mb and 850 mb wind analyses from 0000 GMT 
8 Septenber were used. Comparison of these forecasts shows the forecasts to be 
similar in most areas. However substantial differences in the intensity of 
pressure systems (particularly in the North Atlantic) are evident and these 
differences are enhanced where the initial wind errors extended to the 850 mb 
level . 

The results presented above are typical of other cases and Indicate that 
the forecast model can be sensitive to surface wind data and the effect of such 
data would be enhanced if higher levels were also affected in the analysis. 
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a) Sea level pressure analysis for 0000 GMT 9 September 1978 

b) 1000 mb wind analysis for 0000 GMT 9 September 1978. 

c) 1000 mb wind analysis for 0000 GMT 8 September 1978. 










SEASAT-A IMPACT EXPERIMENTS 
R, Atlas, W. E. Baker, M. Halem and E. Kalnay 


A study Is being conducted to assess the Impact of Seasat-A scatterometer 
(SASS) wind data on GLAS model analyses and forecasts. The SASS wind data 
fields, utilized in this study, were obtained from the Atmospheric Environment 
Service (AES) of Canada, courtesy of Dr. S. Peteherych. The inferred SASS 
winds contain a directional ambiguity with up to four possible directions 
associated with each report. Consequently, this directional ambiguity or 
"aliases" must be removed or accounted for when used in objective analysis. 

The 6LAS system employs a three-pass "iterative" procedure to objectively 
dealias the SASS wind field. Subjective deal lasing and subjective enhancement 
of objectively dealiased winds are also being performed in collaboration with 
Mr. P. Woiceshyn (JPL), Dr. M. Wurtele (UCLA) and Dr. S. Peteherych (AES). 

A preliminary series of assimilation/forecast experiments have been 
performed with the 6LAS analysis/forecast system (previously used for F66E 
experiments). Four assimilation experiments were conducted from 0000 GMT 7 
September 1978 to 0000 GMT 15 September 1978 which differed only with respect 
to the inclusion or exclusion of objectively dealiased SASS wind data or 
satellite temperatures from the available Vertical Temperature Profile Radiometer 
(VTPR) on board the NOAA-4 satellite. A total of sixteen forecasts (four from 
each assimilation) were generated from initial conditions at 0000 GMT 9 September, 
1200 GMT 10 September, 0000 GMT 12 September and 1200 GMT 13 September 1978, 
respectively. 

The SASS data impact is evaluated by means of pressure and wind field 
comparisons, objective measures of forecast accuracy and subjective comparisons 
of prognostic charts. Preliminary results indicate a negligible impact of the 
SASS data in the Northern Hemisphere and a small positive Impact in the Southern 
Hemisphere. However, a larger Impact is obtained in the Southern Hemisphere if 
the VTPR data are excluded from the system indicating some redundancy between 
the two systems. 
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APPLICATIONS OF A GLOBAL SPECTRAL ANALYSIS METHOD 
M. Halem and E. Kalnay 


Current four dimensional analysis (e.g. Rutherford, 1972) combine model 
forecasts with observations of mass and wind fields. In these schemes, the 
model provides an estimate of the fields In regions with data gaps. However, 
the dense coverage afforded by the satellite soundings allows useful estimates 
of not only temperatures but of thermal winds [Petersen and Horn, 1977; Broderick, 
1981], Thus, the uniform coverage afforded by satellite sounders suggests that 
for the first time It may be possible to produce dally global analyses of 
heights and winds without the use of & model. Although such analyses may not 
be optimal, their errors are uncorrelaVed with model errors, which could />rove 
useful In validating results from four dimensional analysis schemes. 

In this work we produce global spectral analyses of height and geostrophic 
winds employing the regularization method of Tyknonov (1963, 1964), In a form 
similar to that of Wahba and Wandelberger (1981). The method Is well suited 
for satellite data but can be used also with other types of data. An optimal 
analysis based on cross-validation, using estimates of observational errors 
Is obtained. For simplicity, the method Is applied first along latitude circles, 
and then 2-d1mens1onal analyses are obtained by expansion In spherical harmonics. 

The global analysis scheme Is shown to be (1) numerically stable and 
convergent In data sparse regions, i.e.. It minimizes errors of Interpolation, 

(11) easy to implement since It avoids the need for solving the partial 
differential equations of mixed type that results In variational analysis 
schemes [Sasaki, I960]. 

Fig. la and lb show the global height field obtained from the NESS and 
6LAS soundings and in Fig. Ic we show the analysis produced from the gridded 
level Ill-b GLAS analysis from OOZ Jan. 5, 1979, all subjected to the same 
spherical harmonic truncation (T-15). In the southern hemisphere, the three 
analyses are Very similar with differences occuring mainly on the Antartic 
continent edge. In the northern hemisphere, the largest differences appear in 
the orientation of the trough over Hudson's Bay, where the GLAS data produced 
a tilt to the NNE while in the NESS system there is a much stronger NE orienta- 
tion. Both systems over Europe are in very good agreement with the FGGE analy- 
sis. 


This method of analysis ought to be particularly applicable to the strato- 
sphere where currently analyses are produced without the benefit of an accurate 
first guess forecast field and where satellites are the major data source. It 
is also a region that does not require high horizontal wave resolution, allowing 
for better fidelity of the harmonic analysis. Figs 2 and 3 compare the 10 mb 
height analysis from NESS soundings with the operational analyses of the United 
Kingdom Meteorological Office (UKMO) and NOAA for two different wavenumber 
regimes. Again, we find these analyses to be In good agreement except near the 
pole. It should be noted that no NESS data were available at 10 mb poleward of 
75°N during the SOP-1. 
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Among the potential applications of this analysis scheme, we mention: 

1) It provides estimates of the state of the atmosphere without the use 
of a numerical model, so that no spurious waves can be Introduced by the model 
during an analysis cycle. This Is particularly important In observational 
studies of atmospheric wave propagation, where free wave« excited In the model 
may contaminate the four-dimensional analysis resultc (Lindzen, Straus and Katz, 
1983). 

2) It provides estimates of geostrophic winds without making use of the 
statistical assumptions commonly Introduced In optimal interpolation, which can 
be grossly In error (Balgovind et al. , 1983). 

3) It can be used to obtain 3-d1mens1ona1 satellite temperature retrievals . 
Current satellite temperature retrievals are one-dimensional, with a column by 
column approach. This does not take advantage of the strong redundance present 
In current radiance observations which are made with a horizontal resolution 

of the order of 10 km, much smaller than even subsynoptic atmospheric scales. 

Clear column radiances should be functions of temperature and moisture, and 
therefore have smooth spatial structure. The present analysis scheme can be 
combined with temperature retrieval methods In the following way: Cloud fil- 

tering schemes are applied on the radiances to produce a first estimate of 
clear column radiances. The variational analysis Is then applied horizontally 
to the clear column radiances. This has two major effects: small scale noise 

can be significantly reduced, and clear column radiances can be estimated even 
for those spots and channels where excessive cloudiness makes cloud filtering 
Impossible. The standard vertical column temperature retrievals are then 
applied on the analysed clear column radiances at a regular or Irregular grid. 

The whole procedure can be performed Iteratively. The net result of this 
method should be to reduce both observational noise and data gaps. 

4) The method can be applied to four»d1mens1ona1 spectral analysis . Halem 
and Kalnay (1983) Indicate that "off-line” global spectral analysis can be used 
as observational Input to four-dimensional spectral analysis. They showe that 
under certain assumptions (I.e., that observational end forecast error statistics 
are horizontally homogeneous, and that In the computation of the error covariances, 
the wave-to-wave Interactions can be neglected), four-dimensional Kalman filtering 
can be drastically simplified. Their analysis shows that with these assumptions 
the formulation of Kalman filtering by 6h11 et al . (1982) reduces to a system 

of scalar equations, one for each wavenumber. 
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APPROXIMATE CALCULATIONS OF GREAT CIRCLE DISTANCES 

R. Hoffman 


For many objective analysis applications it Is necessary to calculate the 
distance between a great many pairs of points on the earth's surface. These 
distances are great circle distances If the approximation that the earth as a 
sphere Is used. (Airlines often follow great circle routes because great 
circles are the snortest path between two points on a sphere along Its surface.) 
It Is the purpose of this note to describe some approximations for calculating 
great circle distances when the two points are nearby. This Is usually the 
case In objective analysis. It Is assumed that the coordinate system Is latitude 
and longitude. Because of the singularity at the poles In this coordinate 
system there Is not one single best solution to the problem. 

A great circle of a sphere Is the Intersection of any plane which passes 
through the center of the sphere with the surface of the sphere. This plane Is 
shown in Fig. 1 ; It Includes the two points. Pi and P2, the great circle and 
the center of the sphere, 0 . If distances are measured In units of the radius, 

R, of the sphere then the great circle distance g Is also the angle P1OP2 
measured In radians. Then by drawing the radius which bisects this angle and 
considering the right triangle which results, we see that 

(1) s1n(g/2) = d/2 

where ^ is the distance (in Cartesian space) between the two points. That Is, 


d2 = (xi - x 2)2 + (yi - yz)2 + (zi - zz)^ 

where (xi, yi, z^) are the Cartesian coordinates of the Hh point. The 
orientation and placement of this coordinate system Is arbitrary. By trans- 
forming latitude and longitude to any convenient Cartesian coordinate system 
It Is tedious but not too hard to show that 

(2) (d/2)2 = s1n2(j» - s1n2x (si 024* - cos2^) 

where 

♦ = (M - *2)72 

X = (Xl - X 2 + 2 nir )/2 
? = (♦! + ♦ 2)/2 
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where X^) are the latitude and longitude of the 1.th point, ji may be zero 
or any Integer here, but In what follows It should be chosen so that [xj jC ir/2. 
Solving (1) and (2) for £ Is relatively expensive as there are five special 
function calls (two sines, a cosine, a square root and an arcsine). 

To get approximations for small £ expressions are found for ^ of the form 


(3) (d/2)2 = r^Cl + a + b + 0(s6)] 


where a = 0(s2), b = 0(s‘*), r = 0(s) and s Is a small parameter. Then by using 
the asymptotic expansions for (1 + x) and s1n“Mx) fs found that 

(4) q = 2rD + {(a 4- r2/3)/2) + {(b + a(rZ - a/2)/2 + 3rV20)/2} + 0(s6)] 


Note that the first term In { } Is 0(s2) and the second Is 0(s**). Three 
cases are considered below for approximations valid away from the pole, near 
the pole and over the whole sphere. The first two cases require two special 
function calls while the third case requires three special function calls; 
this additional call Is the cost of a uniformly valid approximation. 

Case 1. Away from the pole both ^ and x are always small If £ Is small 
so everything but cos ^ ~ a may be expanded. The result Is 


r2 = 1*2 + 0,2x2 


s2 = (j»2 + x2 


a = -(s2 + (2 - a 2 )x 24 , 2 /p 2)/3 


b = (2(s'4 - + (13 - 2a2)s2x2<j,2/r7)/45 


For this case, the leading order behavior, 1,e., neglecting terms of 0(s2) in 
(4), Is the approximation given by Schlatter (1975). In Eq. 4 the more terms 
retained the more accurate the solution as s2 approaches zero but for finite s2 
neglecting terms of 0(s*+) gave the best results In numerical tests. 

Case 2. Near the pole x may be large even though the two points are 
close together. 4> is still small and |<|»| must be close to ir/2. Let n = 

ir/2 - 1^1 so that cos I" = sin n. Then everything may be expanded but sin x = 

M and 

r2 = p2,|,2 + « ii^)n2 


S2 = ^2 + ,^2 
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Case 3., For a qlobally valid approximation neither cos I" = <» or sin ^ » 
P may be expanded, resulting in 


B (1 - + a^P^ 

2 a 2 
S = f 


a = -(1 - p2)4»V3r^ 

b = 2(1 - P^)4>®/45r^ 


In all cases great care must be taken in programming to obtain efficient 
proqi‘ams; all unnecessary arithmetic and logic must be eliminated. Note that 
when r^ = 0 the two points are identical and g = 0. When the approximations 
for cases 1 and 3 neglecting terms of O(s^) were tested against the exact 
solution the relative errors were smaller than 15^ for £ as large as 1, 

Finally, we note that for some applications it may be necessary to cal- 
culate great circle distances between a central point and many other points. 
In such a case one could also make use of the fact that f = <1*1 - <1* is small 
in calculating 


= cos^i^l + 2(cos " 1*1 sin ‘fj) + (sin^i*! - cos^i»i) 

- 4/3 (cos *l’i sin '('i) - 1/3 (sin^+l " cos^+i) <!>'* + ..... 


Reference: 


Schlatter, T. W., 1975: Some experiments with a multivariate statistical 

objective analysis scheme. Mon. Weather Rev>, 103, 246-257. 
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IHPROVEMENT OF MSU TEMPERATURE RETRIEVALS BY USE OF TROPOPAUSE HEIGHTS DERIVED 

FROM TOMS OZONE MEASUREMENTS 

M.-J. Munteanu, E. R. Westwater, and N. C. Grody 


The purpose of the paper is to use the high correlation between TOMS total 
ozone data and tropopause height to improve the MSU temperature retrievals. 

The implementation of the tropopause height information is done in the 
context of stratified regression. 

As discussed by Munteanu (companion paper) tropopause pressure can be 
derived from TOMS measurements of the total amount of ozone. Measurements of 
tropopause height obtained by a completely different technique, namely a ground- 
based VHF radar, led to dramatic improvement of temperature retrievals obtained 
from a ground-based microwave radiometer (Westwater et al ., 1983). It would be 
of large practical importance if similar improvement could be obtained in 
satellite temperature retrievals from a satellite-based instrument which measures 
tropopause height. 

We present here two case studies in Which temperature profiles are derived 
from simulated MSU brightness temperatures, with and without the ozone-derived 
tropopause heigh information. We show that the accuracy and resolution of the 
derived profiles is considerably improved, especially in the region of the 
tropopause. 

Our results are based on high spatial correlation between total ozone 
maxima and minima and the troughs and ridges in the tropopause field obtained by 
M.-J. Munteanu in a pt'evious paper. 

The conclusions of the above mentioned paper are the following ; 

1. The sample correlation was found to be .94 and . 96 , and as a result of 
this high correlation the rms error in tli? prediction of the tropopause 
height from total ozone was found to be 15 mb. In the prediction procedure 
the regression has been used in a manner equivalent to cross-validation. 

2. Correlation between tropopause height and total ozone is the highest among 

all the other correlations with variables such as: brightness temperatures 

from microwave channels, ratio and differences of these quantities, lati- 
tude, mean layer temperatures, surface temperatures or surface pressure. 

3. The high ozone values are associated with lows in 500 mb and 300 mb analysis 
and low ozone values with corresponding highs in the heights analysis. 

4. In the MSU regression retrievals in the low tropopause cases there is a 
large bias and a significant residual error. The tropopause is mislocated by 
100 - 200 mb. 

5. The total rms error in the MSU regression retrievals is rather large in the 
region of the tropopause in comparison with the rest of the retrieval. 
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RETRIEVAL RESULTS 


The results of our stmulatton study will be shown in three complementary 
ways. The first is simply to show typical retrieval results for the three 
coefficient sets (see Fig. 1, 2). All the sets of retrieval coefficients 
used in this study are derived using linear regression. 

The first set is derived from an ensemble of 400 mid-latitudinal profiles 
with 0.5 K noise. We call the regression retrioval obtained global regression. 
The second set was derived from a collection of about 2000 profiles for the 
period February-March-April 1972-1977 for Denver, CO. The regression retrieval 
obtained is called here local regression. The third and most relevant sets of 
retrieval coefficients for this paper were calculated for ensembles containing 
profiles stratified by tropopause pressure within 40 mb thick class intervals. 
The data base for the tropopause stratification was all profiles 1972-1977, 
Denver, CO. 


CONCLUSIONS 

As is seen from the Figs. (1, 2) a substantial improvement in accuracy occurs 
around the tropopause region when stratified coefficients are used. In all of 
our cases, no degradation of retrievals occured when tropopause stratified 
coefficients were used, and in a substantial fraction cases, a dramatic improve- 
ment in accuracy and resolution was apparent. 

Sometimes even in the lower tropopause region the accuracy of the retrieval 
is improved. 

The general improvement in accuracy obtained by using ozone data is quite 
apparent in the statisical sense. In Table 1, we present at mandatory pressure 
levels the rms errors of the three sets of retrievals for the April 15, 1979 
and the April 25, 1979 cases. The performance of the two climatological 
retrievals is bettered at almost all levels by the tropopause stratified 
retrievals, with improvements ranging from 1 to 2 K rms. 

As a final illustration of the retrieval improvement offered by the utili- 
zation of ozone data, we present two-dimensional fields of tropopause height 
as derived from (a) radiosondes (ground truth), (b) local climatological retrie- 
vals, (c) global climatological retrievals, and (d) tropopause stratified 
retrievals. It is apparent that both climatological retrievals result in 
tropopause height fields that bear little resemblance to the actual ones. 

However, the thermal fields derived from the combined MSU and ozone contain 
several, but not all, of the features of the original field. These results 
strongly suggest that tropopause height information is, at best, only weakly 
contained in the MSU radiances. 

As a general conclusion the results presented strongly suggest that a 
substantial improvement in, at least, MSU thermal retrieval accuracy can be 
obtained by utilizing TOMS measurements of total ozone content and its subse- 
quent relationship to tropopause height. These results are totally consistent 
with earlier efforts incorporating radar tropopause height measurements with 
passive soundings (Westwater et al ., 1983). 
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STUDY OF TROPOPAUSE HEIGHT ESTIMATE FROM TOMS TOTAL OZONE DATA FROM 
NIMBUS-7 AND FROM THE MICROWAVE REGRESSION TEMPERATURE RETRIEVAL OF SIMULATED 

BRIGHTNESS TEMPERATURES 

M-J. Munteanu 


A well known characteristic of satellite soundings is the fact that they 
exhibit very large errors in the region of the tropopause due to poor vertical 
resolution (N. Phillips et a] . , 1981, Schlatter, 1973, NESS group, 1982, Susskind 
et al . , M. Halem et al .,T98T). 

This study pursues the idea of investigating the potential use of TOMS 
total ozone from NIMBUS-7 satellite in order to improve the satellite retrievals 
in the region of tropopause. 

We have concentrated our research in the following four areas; 

(i) Determination of the tropopause using TOMS data. 

(ii) Impact of tropopause errors in microwave MSU retrievals. 

(iii) Impact of tropopause errors in the future AMSU retrievals. 

(iv) Use of air-mass to stratify regression retrievals. 

( i ) Determination of the tropopause height using TOMS data . 

The first step has been the investigation of whether an accurate estimation 
of the tropopause height can be obtained using total ozone remotely sensed from 
NIMBUS-7 satellite. For this purpose we have studied the correlation between 
TOMS total ozone and tropopause height from radiosonde data over Europe for 
two areas of interest. 

I Northern Europe; 10°E - 50°E and 50°N-65°N Date: April 15, 1979 

II Western Europe: 10°W-20‘>E and 30°N-65®N Date: April 25, 1979 

The data consisted of NIMBUS-7 TOMS total ozone data colocated in time and space 
with radiosonde measurements of temperature and humidity and NMC geopotentlal 
heights analysis. The conclusions from this study are: 

1, The sample correlation between tropopause height and total ozone data were 

found to be .94 and 0.96 which corresponds to an rms error of about 15 mb. 

For Comparison the prediction of the tropopause height from radar measure- 
ments has an rms error of 30 mb. In the prediction procedure the regression 
has been used on the manner equivalent to cross-validation. 

2. Correlations between tropopause height and total ozone is the highest com- 
pared with all other correlations with variables such as: brightness temp- 

eratures from microwave channels, ratio of differences of these quantities, 
latitude, mean layer temperatures, surface temperatures or surface pressure. 


3. Experiments correlatinq tropopause data with mean layer derived directly 
from radiosonde data showed that even predictions of tropopause height from 
exact mean layer temperatures were inferior than those derived from total 
ozone. 

4. The folding of the tropopause should be taken into account because it is 
associated With the highest gradient on the ozone field. We should give 
different weights to the prediction in those regions where the folding of 
the tropopause is likely to occur. 

5. The high ozone values are associated with lows in the 500 mb analysis and 
the low ozone values with corresponding highs in the heights analysis (see 
Fig. D'J). 


( i i ) Impact of tropopause errors in the MSU retrievals . 

To study the impact of the lack of resolution at the tropopause levels we 
have performed simulation of microwave radiation measured with three MSU chan- 
nels existing on TIROS-N and then performed local regression retrievals with 
Westwater's formulas for April 15, 1979 and global regression retrievals with 
Grody's formulas for April 25, 1979. 

The conclusions are: 

1. In the low tropopause cases there is a large bias and large residual error 
(Figure [2]). 

2. The low tropopause cases produce very poor temperature retrievals. The over- 
all rms is large in the region of the tropopause (Figure [2], Figure [6]). 

3. The tropopause is mislocated by 100-200 mb in the regression microwave 
retrievals in the extreme cases. In one example the tropopause is esti- 
mated at 300 mb and it is actually observed at 490 mb. In another case 

the tropopause is observed at 450 mb, 419 mb, 412 mb, 415 mb, and estimated 
in all these cases at 300 mb (Figure [2]). 

Two dimensional fields of: ozone, tropopause height, predicted tropopause 

height from total ozone, tropopause height calculated from climatology regression 
retrieval are presented in Figure [3] and Figure [4] for the 15th and 25th 
respectively. 

4. As it is observed the MSU regression temperature retrieval damps almost 
entirely the trough in the two dimensional field of the tropopause height 
for April 15 and in general does not contain all of the highs and lows for 
April 25. 


( i i i ) Impact of tropopause errors in the future AMSU retrievals 


Experiments involving microwave simulation of radiosonde data with 5 of 
the AMSU channels and a comparison with MSU simulation have been performed. 
The conclusions are the following: 

1. The location of the tropopause is not improved with respect to the MSU in 
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the low tropopause cases (see Figure [2]). 

2. The retrievals based on the 5 AMSU channels generally perform better in the 
average tropopause height cases. 

3. The general conclusion of these experiments could be described by saying 
that the resolution in the region of the tropopause is not improved by the 
additional AMSU channels, especially in the cases of low tropopause. 


(iv) Use of air-mass to stratify regression temperature retrievals 

The correlation between the tropopause height and the difference between 
the brightness temperatures of the two MSU channels situated in the region of 
the tropopause is greater than with any other combination of individual chan- 
nels. This leads to the natural idea of using the difference between the two 
MSU channels as a stratifier for regression temperature retrieval formulas. 

The main conclusion is: 

1. Experiments using stratified temperature retrieval formulas indicated that 
they perform better in the region of the tropopause (see Fig. [5]). 
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THREE DIMENSIONAL INVERSION OF SATELLITE OBSERVED BRIGHTNESS TEMPERATURES 

R. N. Hoffman 


We propose here to obtain fully three dimensional fields of temperature 
from observed satellite brightness temperatures by using some of the variational 
techniques championed by Wahba and collaborators (Wahba and Wendelberger, 1980; 
and references therein). We recognize from the start that radiances or the 
equivalent brightness tt/iperatures are not tropospheric temperatures but only 
proxies for these temperatures and that our results will depend upon what a 
priori Information the Inversion process utilizes. ThTs being said It folTows 
that the best combination of all a priori Information should be used. The 
variational procedure allows In a natural way a variety of different Informa- 
tion to contribute to the Inversion process. In calculating three dimensional 
retrievals we can fully use our a priori knowledge that horizontally the tem- 
perature field Is smooth to obtain maximum vertical resolution, filter cloud 
effects, minimize the effects of Instrument noise and obtain horizontally 
coherent retrievals. 

It has often been remarked that satellite observed radiances contain more 
Information about temperature gradients than Information about temperature. 
Because it acts as a spatial filter, there are several reasons to expect that 
the var1at1onf«l analysis method (VAM) will result In better and more useful re- 
trievals. First, the retrieved temperature field will be horizontally coherent. 
The effect of random Instrument noise and small scale meteorological features 
which are not really resolvable by the Instrument will not be present In the 
retrievals. Second, by combining Information from many spots (that Is obser- 
vation points or satellite footprints) the maximum amount of Information about 
the vertical structure of the atmosphere may be extracted. This Includes 
cloud structure; the VAM has the potential to act as a very powerful cloud 
filter retrieving cloud amounts, heights and other properties (such as emlssl- 
vlty). Third, the VAM Is very flexible In allowing additional sources of 
Information such as a 6 h forecast or the conventional radiosonde data to 
influence the retrieved temperatures. In addition the VAM may be fully Inte- 
grated with a NWP model to provide what might be termed four dimensional retrie- 
vals. Fourth, the above points also apply to retrieved surface quantities 
such as temperature and emissivity. 

Notwithstanding the above observations, every retrieval system attempted 
to date has been essentially one dimensional. (By one dimensional we mean that 
the clear column radiances for each spot are processed independently of each 
other spot. Of course, for IR retrievals, neighboring spots are used to remove 
cloud effects and thus reduce the 'raw' radiances to clear column radiances). 

One might suppose that If the one dimensional problem is difficult the three 
dimensional problem will be more difficult and should not be attempted until 
the one dimensional problem Is solved. It Is certainly true that the three 
dimensional problem Is very complicated and difficult computationally. However 
the three dimensional problem may be better posed In the sense of uniqueness 
and stability. 

The basic problem of three dimensional retrievals, common to many Inverse 
problems (geophysical, medical, etc.). Is that there Is a finite number of 
observations of some parameters which depend on the continuous distribution of 
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some properties of the system under observation. In the present case we observe 
brightness temperatures 




at a number of viewing geometries (Rn) in a set of channels i^i), and these 
brightness temperatures depend on among other things, mixing ratios of the 
absorbing gases and the temperature all along the path followed by the radia- 
tion through the atmosphere, The result is that the problem of retrieving the 
system properties from the observations is grossly underdetermined and to a 
large extent the solution depends on what other information is added by the 
investigator (Backus and Gilbert, 1967), In many cases it is possible to 
avoid this problem simply by limiting the resolution of the solution, that is 
to use a relatively coarse grid to represent the system. This presupposes 
that the system properties vary smoothly. This assumption is in fact addi- 
tional information. In the atmosphere temperature is very smooth, the mixing 
ratio of water is smooth but less so and cloud amount is sometimes not smooth 
at all on the length scale of the observing cell. 

The difference between the length scales for temperature and cloud amount 
is the crux of the three dimensional retrieval problem. In general we can 
classify a system property as either ragged (resf'!» smooth) if its length scale 
is the order of or smaller than (resp. greater than) the observing cell size of 
the instrument. For properties which are ragged, like cloud amounts, we may be 
able to infer some sort of average value over each particular observing cell 
but it is impossible to retrieve point values. Smooth properties on the other 
hand should be retrievable at points either within or between the observing 
cells. It is important to note that while cloud amounts are ragged, cloud 
properties (such as emissivity) are smooth if we allow their definition where 
there are no clouds (Fig. 1). This is the basis of the multiple fields of 
view cloud filtering methods (Chahlne, 1977; and references therein). The 
properties of the earth's surface are also in general ragged because the length 
scale for soil types and land use is relatively small. 

The intuition behind our strategy is therefore the following: We expect 

to use one channel for each ragged property needed for the forward problem. 

Then the remaining channels will be pooled geographically by the smoothness 
constraints acting on the smooth parameters to yield optimal estimates. Of 
course the better we know the smooth properties the better we can estimate the 
ragged properties and vice versa. Within the variational formulation the 
retrieval of all quantities occurs simultaneously and this beneficial interac- 
tion between knowledge of the various variables will occur naturally. 

To generalize the variational problem of Wahba to take all relevant para- 
meters into account let us symbolically write the solution to the forward 
problem as 


^c(^i* ^n ^ “ ^i (pi(^n)* P2(^n)***Pr^^n)* ®l(^n)***®s(^n)) 


where Pi(Rn) represents the value(s) of the j^tn rough variable entering the 
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forward problem evaluated for viewing geometry R^; similarly aj(Rn) represents 
the value(s) of the jth smooth variable entering the forward problem. Some of 
these variables are 2-dimensional such as surface emissivity, while others are 
3-dimensional like air temperature. In general some of the pj will be cloud 
amounts for each of several cloud types and some of the aj will be the corre- 
sponding cloud top pressures, emissivites, etc. Now we wish to minimize 


(1) S s* y Win (Tb(vi, Rn) - Tc(vi, Rn))^ 

1 »n 

s 

j=l J 

where 

m! 

'^m(®) “I ^ aj !«2 !« 3 ! ^ /// ^ 3x**^3y“^3p“^ dxdydp 

aj^+a2+a3=m 


During the minimization m and X are held fixed. We must still specify f, the 
solution of the forward problem, the spatial representation of 0 and the 
remaining free parameters (m, X, w^p) to complete the problem. Note that 
a different way of measuring the smoothness of each of o-variables is possible 
and that there is no penalty if the p-variables are not smooth. Also there 
is a X weight for each o-variable to balance the tradeoff between the various 
constraints. One could, and probably will, append to the end of (1) a variety 
of other constraints, such as climatological information, 'first-guess' infor- 
mation, radiosonde data, etc. Although the mathematical machinery developed 
by Wahba ,^pplies directly only to linear functionals, the general method has 
already been applied to large nonlinear problems (Hoffman, 1982). 
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Fig. 1. Cloud parameter (p) versus position (r). Cloud 
amount (n) is ragged and Boolean. Cloud emissi-^ 
vity (e) is smooth within clouds and may be 
defined smoothly between clouds. 
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SIMULATION COMPARISON STUDY OF THE AMTS AND HIRS2 SOUNDERS 
J» Susskindj D. Reuter and A. Dalcher 


The Advanced Moisture and Temperature Sounders (AMTS) is a proposed high 
spectral resolution passive infra-red sounder designed to produce higher vertical 
resolution soundings of temperature and humidity than currently available 
[Kaplan et al .. 1977]. In addition, the channels are chosen so as to essentially 
eliminate any dependence of the signals on extraneous variables, i.e., H 2 O and 
O 3 distribution in the case of temperature sounding channels. 

The simulation study is part of a joint study between NASA, NOAA, and the 
University of Denver to assess the sounding capabilities of the AMTS as compared 
to the HIRS2 instrument currently flying on the TIROS-NOAA series operational 
satellites. Radiances for both instruments were simulated by the University of 
Denver for a set of profiles containing information of the ground temperature, 
atmospheric temperature-humidity-ozone profiles, reflected solar radiation off 
the ground, and satellite zenith angle of observation. In one portion of the 
test, radiances were simulated for a large sample of cloud free conditions. In 
a second portion, more detailed simulations were made involving radiance fields 
computed as a function of multi-layer cloud fields, thermal gradients, and the 
scan pattern of the instrument. Only the clear test is considered in this 
write-up. The cloudy portion is discussed in Dennis Reuter's section in this 
report. 

In the clear test, radiances were simulated by the University of Denver 
for 1600 profiles, selected from 3 week periods in each the winter and summer. 

The radiances for both instruments, together with colocated radiosonde informa- 
tion, were sent to both NOAA and NASA. Each group analyzed data for each 
instrument according to their preferred method of analysis, regression analysis 
in the case of NOAA, and a physical inversion method by NASA, The test then 
compared both methods and instruments. The colocated data was to be used for 
any purpose i.e., generation of regression coefficients for NOAA and generation 
of empirical tuning coefficients to minimize systematic differences between 
radiances calculated by the University of Denver and by NASA for the same pro- 
files. The actual test was based on radiances simulated from radiosonde data 
in the subsequent 2 week period. These profiles were not made known to the 
participants. Here, we briefly show the GLAS method of analysis and indicate 
results based on the original sample of 1600 cases, in which half the radio- 
sonde information was withheld in the generation of coefficients and that 
sample used as the test data set. 

The method of analysis was very similar to that used by GLAS in analysis 
of TIROS-N HIRS2/MSU data [Susskind et al., 1982]. The procedures used in the 
test differ from those used in analysis of TIROS-N for two reasons: 

1) A 6 hour forecast guess of temperature and humidity is available in analy- 
sis of TIROS-N data but is not available in the test. 

2) The noise levels in the test are realistic assessments of instrumental 
noise but do not include the effects of scene noise. 
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As a result of these differences, the first guess used in the analysis was 
based on a regression relationship between observed brightness temperatures and 
radiosonde temperature profiles. Also, the form of iterative relaxation equa- 
tion was modified so as to decrease the smoothing applied to the solution and 
a correction was applied to the climatological humidity as a first guess pro- 
file. 


In order to generate an initial guess temperature profile based on regres- 
sion relationships between observed brightness temperatures and atmospheric 
temperature profiles, it is desirable to remove to first order the angle depen- 
dence of the satellite brightness temperatures. To do this, observed satellite 
brightness temperatures at angle 6, Tb( 6), are corrected to Tb( 0), their estimate 
values if the observations were at nadir. 


TABLE 1 

CHANNELS USED IN DIFFERENT STEPS 



HIRS 

AMTS 

ANGLE CORRECTION 

1-8, 13-16 

3-10, 20-26 

REGRESSION GUESS 

1-8, 13-16 

3-10, 20-26 

TUNING 

1-7, 13-16 

1-10, 20-26 

GROUND TEMPERATURE 

18, 19 

26-28 

HUMIDITY CORRECTION 

8 

11 

TEMPERATURE PROFILE 

1-4, 13-15 

4-10, 20-24 
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TABLE 2 

ASSIGNED PRESSURES FOR CHANNELS USED IN THE TEMPERATURE RELAXATION SCHEME 


AMTS HIRS 


CHANNEL 

Pi 

PL“Pu 

CHANNEL 

H 

Pl-Pu 

9 

3 

- 

1 

10 

- 

10 

15 

- 

2 

- 

30-90 

8 

- 

30-90 

3 

- 

90-200 

7 

- 

30-80 

4 

- 

200-380 

6 

- 

50-150 

15 

- 

380-625 

5 

- 

100-220 

14 

- 

625-875 

4 

- 

200-400 

13 

- 

875-1000 

20 

- 

300-500 




21 

- 

400-600 




22 

- 

600-775 




23 

- 

775-1000 




24 

* 

925-1000 





We use the equation 

SGC 0-1 

i(0)=T|(6)+A T§(0) (1) 

sec 50-1 


The superscript S means simulated data. A is 15 x 15 Tor AMTS and 12 x 12 for 
HIRS2 as shown in table 1. A separate A matrix was constructed for winter and 
summer. The coefficients of A were determined from simulated radiances for the 
800 profiles at a 50° zenith angle and at nadir. 

Once the Goldman data were corrected to zero zenith angle by equation (1) 
regression equations of the form 


T(P) = T(P) t B[T«(0) - Tb(0)] 


( 2 ) 
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were constructed for each of eight zones, (land/water), (tropics/extratropics), 
and (summer/winter). Tb( 0) is the vector of brightness temperatures for the Nr 
channels used for regression constructed by angle correcting the data provided 
by the University of Denver for the particular profile to 0“ according to 
equation (1). Tb( 0) is the mean vector of all the resulting Tb( 0) in the data 
set. Nr = 15 for AHTS and 12 for HIRS with the channels shown in table 1. All 
humidity sounding channels and most window channels have been excluded from the 
regression for both instruments so as not to use humidity as a predictor for 
temperature. 

Possible systematic computational differences are minimized by comparing_a 
set of "observed" brightness temperatures To, with those simulated by GLAS, T0, 
and finding C and D which best fit the relation 


= C Tg + D (3) 


where C is an Nj x Ny matrix for the Nj channels tuned and D is a Ny x 1 vector. 
Once C and D are obtained from analysis of a dependent set, the observed 
brightness temperatures Tg in the independent set are modified according to 


Tb = C Tb + D (4) 


so as to best match the brightness temperatures we would compute under the same 
conditions. 

Retrievals are performed using Tb as data rather than Tb. We now compare 
the modified observed brightness temperatures for the temperature sounding 

channels Tg i, with the computed brightness temperatures from the Nth guess, 

Td j , and modify the guess to produce an N+lth iterative temperature profile. 
Twelve temperature sounding channels are used for the AMTS and seven are used 
for HIRS2 as shown in Table 1. 

The relaxation equations used are almost identical to those described in 
Susskind et al . [1982]. Temperatures at pressures between 30 and 1000 mb are 
treated differently than temperatures at pressures lower than 30 mb. The basis 
of the relaxation method lies in the approximation that small a constant shift 
in the temperature profile will produce an almost identical change in the 
brightness temperature computed for a sounding channel. Moreover, if the shift 
is applied only in the region of the atmosphere where the radiance of the 
channel is most sensitive to atmospheric temperature, a similar change in 
computed brightness temperatures will occur. 

For channels sensitive to temperatures at pressures greater to 30 mb, we 
assign an atmospheric layer, PiL To Piu, which is considered to be the layer in 
which the radiances of the channel are most sensitive to atmospheric temperature 
change. P^l and P|u» representing the lower and upper pressure boundaries for 
sounding channel i, are shown in Table 2 for the sounding channels of AMTS and 
HIRS used in the analysis. Channels primarily sensitive to temperature changes 
above 30 mb are treated as representative of temperature changes at a specific 
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pressure, rather than in a layer. This pressure, Pi* is indicated in Table 2 
for the appropriate channels. 


For pressures greater than 30 mb, we write the relaxation equation 




+ TB.i 



where T^ is the average temperature of the Nth iterative guess temperature 
profile in the atmospheric layer corresponding to channel i 


( 5 ) 



PiU 

/ tI^CP) dtnP /tn 

PiL 


PjU 

PiL 


(fi) 




and ij ls 
atmospheric 
at pressure 
solution to 


the new estimate of the layer mean temperature in the appropriate 
layer. In order to determine a iterative temperature profile 

Pji from estimates of layer mean temperatures, we constrain the 
be given by 


(P|() ' T®(P|<) tZ a5*Sj(P|<) (7) 

where T^(P|^) is the initial guess, Fj(P|^) are empirical orthogonal functions of 
temperature, given by the eigenvectors, with largest eigenvalues, of the cova- 
riance matrix of a set of global radiosonde profiles, sampled at the 52 pres- 
sure levels between 1000 and 30 mb, which are a subset of the 64 pressure 
levels used in the calculation, and A^+1 are the iterative coefficients 
which, together with the initial guess, completely determine the solution. 

The coefficients, are solved for according 


AN+1 = [F*f + oH]-l F'[T^+1 - T^] (8) 


where F represents the matrix of layer averaged empirical orthogonal functions 


_ PiU , ^ PiU 

= / Fj(P) dtnP /An , (9) 

PiL PiL 

F' is the transpose of F, H is a diagonal matrix with Hji being the inverse of 
the fraction of the total variance arisj^ng from eigenvector j, and a is a 
constant. The term oH is added to F'T in order to stabilize the solution, as 
in Susskind etal. [1982]. 
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For AMTS, 10 mean layer temperatures, shown In Table 2, are used to estimate 
coefficients of 9 empirical orthogonal function, with o » 1 x 10“3, For HIRS2, 

6 mean layer temperatures are used to estimate coefficients of 5 empirical 
orthogonal functions with o = 5 x 10"^. 

Figs. 1 and 2 show the RMS errors of the regression generated guess are 
the retrieved temperature profiles obtained from both the AMTS and HIRS data 
for the Independent sample. Separate statistics are shown for the winter and 
summer cases. Only levels which were interpolated from radiosonde reports are 
Included In the statistics. Under clear conditions, the RMS errors of the 
physical retrievals for both AMTS and HIRS are .1® - .2® better than their 
regression guess beneath about 300 mb. Above 300 mb, the physical retrieval 
Is of comparable quality to the guess for the AMTS and somewhat worse for the 
H1RS2. The AMTS results show significant Improvement over the HIRS2 results, 
relatively Independent of the method used. Improvements are the order of .3® 
below 300 mb and .7° above 300 mb. Improvements In the lowest atmospheric 
level, from 880-1000 mb, are the order of .5®. 

Figure 3 shows similar statistics for the single zone land-mid-latitude 
winter, characterized by large Inversions In the lower troposphere. Under 
these conditions, the physical retrieval method shows larger Improvements over 
regression In the lower troposphere and likewise, the Improvements In the 
retrieval actuary obtained frotn the higher vertical resolution of AMTS over 
HIRS are more significant. 
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Figure 1. Results of analysis by GLAS of 400 sets of radiances from the clear 
winter dependent set using information from the other 400 profiles. 
The errors plotted are of the layer-mean temperatures with layers 
shown in the figure. In each case, the error for the layer is 
plotted at the lower boundary of the layer. 
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Figure 2. As In Figure 1 but for the clear summer profiles. 
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SIMULATION STUDIES OF THE HIRS2/MSU AND AMTS/MSU SATELLITE SOUNDING UNITS: 

CLOUDY CONDITIONS 

D. Reuter, J. Sussklnd and A. Daleher 


Currently* joint NOAA/NESS - NASA simulation study is under way to I 

compare the accuracies of atmospheric temperature profiles retrievals as pro- | 

duced from the currently operational HIRS/MSU sounding system and the proposed j 

AMTS/MSU sounding system. This study is proceeding in two phases. In the t 

first phase, which has already been completed, synthetic satellite radiances j 

are calculated from realistic atmospheric temperature, water vapor, and ozone ! 

profiles, with the effects of surface temperatures, satellite viewing angle. If 

and reflected solar radiation also taken into account. These simulated radiances | 

are then used as input data for the statistically based temperature retrieval l| 

algorithm at NOAA/NESS and for the physically based retrieval system at GLAS. | 

This first phase is described elsewhere in this report (see Susskind et al .) | 

and will not be discussed here. |i 

I 

In the second phase of the test, the simulated radiances are contaminated I 

by the effects of realistic cloud fields. In order to prepare for this phase j| 

of the test and, at the same time, to develop and test various methods for use I 

With real data, we at GLAS have carried out simulation studies employing cloud 
contaminated radiances generated here. It is the purpose of this note to 1 

describe these studies. | 

For the purpose of these simulation studies, and in fact for the joint 
test, the clear column radiances are calculated according to a truncated version | 

of the radiative transfer equation: | 


tnF dTi(P,e) 

Ri = BiCTs]Tis(Ps.0) + / BiCT(P)] (1) 

tnPs dinP 

where R^ is the frequency averaged rad farce in sounding channel i, B|[T] is the 
Planck function evaluated at frequency Vj and temperature T, t^(P,0) is the 
channel averaged transmittance from pressure P to the satellite pressure P at 
angle e, and the subscript S stands for surface. This form of the equation 
neglects effects of reflected solar and downwelling thermal radiation but 
preserves most of the dependence of the radiances on atmospheric and surface 
properties and viewing angle. 

In the presence of clouds, which are essentially opaque to infrared 

radiation, the radiance observed by the satellite, R^ , becciTies a function of 
the cloud top pressures Pj, and fractional amounts oj: 

R( = (1 - Z«j) RP-" ♦ I Rjf-® (2) 

J J 

where R^^^ is the radiance which would be observed in the absence of clouds 
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and is the radiance which would be observed if the field of view of the 

sounder were obscured by a cloud with a cloud top pressure of Pj. Equation 2 

Is valid for any value of the emissivity of the clouds, however for the purposes 

of test the emissivities were taken to be.l (i.e. the clouds are taken to 

be perfect black body emitters). Note that the cloud fractions aj are 
effective cloud fractions in the sense that if a high cloud obscures a Tower 
cloud, then only that part of the lower cloud which is not obscured will 
contribute to the radiance. 

If the set of radiances is used to determine the surface and atmospheric 
properties the results will be grossly in error. Therefore some method must be 
used to account fcr effects of clouds. At GLAS we use a method whereby the 
radiance which would be observed in the absence of clouds is reconstructed from 
observed radiances in several nearby fields of view. 

An iterative estimate of the clear column radiance, may be 

obtained from two fields of view as: 


rCLR(N) . + „(N) ,3) 

where R^ j is the observed radiance for channel i in field of view j and is 
a channel'^ independent cloud filtering parameter which is equal to 01/(01-02) 
where 01 and 02 are the effective cloud fractions in field of view 1 and 2 
respectively. 

.In the GLAS physical inversion method we always have an estimate of 
I^LLR^N) iterative atmospheric and surface properties. Therefore 

on any iteration may be estimated as: 


n(N) 


R,«'-C(N) . 




Ri ,1 " Ri ,2 


(4) 


Since it is not expected that, at the beginning of the iterative temperature 
retrieval procedure, the N^h solution profile is in fact the correct atmospheric 
temperature profile, this method will tend to give erroneous results. To 
remove, at least to first order, any biases in rvALCin)^ microwave channel 2, 
which is not affected by the. presence of clouds, is used to obtain a corrected 
clear column radiance 

Rf (N) = B,[t5'''-C(N) + (5) 


where 1s the black body temperature for Is the 

iterative brightness temperature of microwave channel 2 and Tu2 is the 
observed microwave channel 2 brightness temperature. Brightness temperatures 
are used^because the radiances are entceinely frequency dependent. /Ibis value 
for Rp''^rmay be substituted for r^alunj equation 4 to find n'^' which 
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may then be used in equation 3 to obtain the reconstructed clear column radiance 
Since n is channel independent all infrared radiances may be 
corrected by the same n. 

This procedure is quite adequate for accounting for one cloud level and 
therefore one degree of vertical inhomogeneity, however, in the simulation test 
described above multiple cloud levels are to be used and therefore multiple 
degress of inhomogeneity must be accounted for. We have therefore attempted to 
use three fields of view to account for at least two degrees. of inhomegeneity. 

To determine the two cloud filtering parameters and channels 
7 and 6, are used for the HIRS instrument while channels 1 and 2 are used for 
the AMTS. The observations are sorted into warm, intermediate and cold fields 
of view (Fi, F2, F3) according to the brightness temperatures in the window 
channels 8 for HIRS and 11 for AMTS. For three fields of view equation ( 3 ) may 
be expanded to: 


R^CH) , t n{N) (R,| 1 . R,^3) + - R,,2). 


( 6 ) 


Note that is being used to describe most of the inhomegeniety since, 
by the field of view sorting given above (R^^ - Ri,3) >. (Ri,l - Ri,2)* 

We will describe the determination of and by referring 
to the HIRS instrument. If only channel 7 were used to determine m we would 
use equation 4 to obtain: 


n 



r^CCCN) . 


•^ 7,1 


R7,1 - R7,3 


( 7 ) 


A similar equation would be used to obtain m if only channel 6 were used. 
We shall take n|*'*^ to be a weighted average of and 01^5: 

n{*^^ * 4,7 (*^7,1 - *^7,3)^ ^ 4!e (*^6,2 ■ *^6,3)^ 

(^7,1 " *^7,3)^ (*^6,2 - *^6,3)^ + ® 


where 6 is a damping parameter on the order of the microwave noise which is 
analogous to the ridge parameter in ridge regression. This equation is the 
least squares solution for if equation 6 is used to obtain the best 

estimate of based on channels 6 and 7 . In other words and n|^5 
are weighted by the amount of signal in channels 7 and 6 respectively. ’ 
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To determine ng a similar procedure is followed except that part of the 
inhomogeniety described by ni is removed from the signal, for example: 


4 ^)- 


- R7^i) - ^R;^l - R7^3) 

■ II— ■ II.I- ■ I ■ I Wi IIM H |IM» 111— ■ vy 

A/ A/ 

(R7,1 - R7,2) 

,2 


^ 4!’ ("7,1 - "7.2)^ * "M (^6,1 - "e.a)’ 

A4 fW AW />/ 

(R7,i - R7.2)^ + (Re.i - Re, 2)^ + « 




( 10 ) 


If the denominator of equation 10 is found to be on the order of the 
expected combined infrared and microwave instrumental noise n£ is set to 
zero. Similarly if the denominator of equation 8 is on the order. of the expected 
noise is set to zero. .^If is found to be less then zero is set 
to zero, in which case is set to zero since nj”’ should contain most 

A# A/ A> A4 

of the cloud information. If either [(17.1 - 17^3) + (Tg^i - T6,3)]/2 _< 1° or 

/>j r>J A/ rw 

[(^7,1 “ T7 2) + (^6,1 " T6^2)3/2 < 1° it Is assumed that field of view 1 is 
clear and therefore ni and 02 are both set to zero. In this case a further 
check is made to determine if the scene is completely overcast. If - 

^7,1) 2. scene is assumed to. have complete cloud cover and no retrievals 

are done. Finally, if [(n^ + _> 4 the scene is assumed to be too 

cloudy to perform a retrieval. 

The simulation retrieval studies to be described below were carried out 
using cloudy radiances simulated at 6LAS. Retrievals were carried out using 
radiances simulated from radiosonde atmospheric temperature profiles at various 
satellite viewing angles. The surface temperature was taken to be randomly 
distributed about tli* surface air temperature with a standard deviation of three 
degrees. Equation 2 was used to construct cloud radiances for two cloud layers 
randomly distributed between 300 and 700 mb. For each radiosonde three fields 
of view were simulated. The two cloud top heights were taken to be the same 
for all three fields of view, however the cloud fractions for the two cloud 
heights varied randomly for the three fields^ of view._ The avera^ge total cloud 
fraction for the three fields of view were 01 = 0.2, <*2 “ ^.5, a3 = 0.75. 

Thus the cloudy radiance simulations produced, on the average, very cloud 
contaminated radiances. These simulated cloudy radiances were contaminated by 
reasonable instrumental noise and then used as observations for the GLAS physical 
retrieval scheme. The retrieval system then precedes by the steps: 

1) Calculate clear column radiances based on a climatological first guess 

2) Use equation 7, 8 and 10 to calculate nj(^) and 02^^^ 

3) Reconstruct clear column radiances using equation 6 


71 


4) Obtain a ground temperature 

5) Use the new ground temperture to recalulate the clear column radiances 
and use steps 2 and 3 to reconstruct clear column radiances 

6) Perform a regression to obtain an improved temperature profile 

7) Recalculate the ground temperature 

8) Reconstruct clear column radiances 

9) Perform a physical retrieval step and check differences between calcu- 
lated clear column brighness temperatures and reconstructed clear 
column brightness temperatures for convergence 

10) Obtain a new ground temperature 

11) Iterate steps 8, 9, 10 until convergence is achieved (step 10 is only 
done 3 times no matter how may iteration are done for steps 8 and 9) 

12) Check rejection criteria and flag the retrieval if it should be rejected 
(i.e., if the RMS difference between the computed and reconstructed 
brightness temperatures is greater than 1° or if the calculated and 
observed values of MW channel 2 differ by more then, 1°, or 1.5" if both 
ni and 03 = 0) 

Figure 1 shows a comparison of the results of the retrieval for the HIRS2 and 
AMTS sounders. Also shown in this figure are the results obtained from the 
regression analysis performed in step 6. The regression matrices were obtained 
from clear column brightness temperatures obtained from calculations using 
radiosondes not used in this study and thus these profiles represent an 
independent set of regression results. As may be seen from this figure the 
AMTS improves the HIRS results throughout the atmosphere. Furthermore the 
physical retrieval method shows much improvement over the regression results. 

In this connection it should be noted that attempts to improve the regression 
results by iterating the regression (i.e., by substituting step 6 for step 9) 
gave rise to some improvement of the regression results, but they were still 
inferior to the results from the physical retrieval. 

Figure 2 compares the results obtained from these cloudy radiances with 
results obtained when the same radiosondes were used to simulate radiances with 
no clouds. As may be seen from this figure the degradation of the results in 
the presence of heavy clouds is not very severe. 

Finally, to indicate the relevance of these results to the joint NOAA/NESS- 
NASA Simulation Test, Figure 3 compares the results of retrievals performed in 
the clear column portion of this test using radiances generated at the University 
of Denver and the results for the same profiles using clear column radiances 
simulated by the model used in the 6LAS simulations. As may be seen from this 
figure the results obtained using radiances based on the GLAS simulation model 
are quite similar to the results obtained from radiances based on the model to 
be used in the actual test. 


As a final note it might be asked what relation, if any, this simulation 


study bears to what might be expected from an operational sounder sensing the 
true atmosphere. Since the GLAS physical inversion method relies quite heavily 
on the ability to calculate the radiances which would be observed at the satel- 
lite given the true atmospheric and surface conditions, the errors in this 
calculation will very probably be the limiting factor in determining the appli- 
cability of these results to real data sounding. However, it has been shown 
(Susskind et al ., 1982) that the model used in this study quite accurately 
reproduces satellite observation when the relevant parameters are known. 
Therefore, we feel that we may transfer results obtained from this work to the 
problem of actual meteorological parameter sensing. In particular the current 
analysis system used at GLAS for the HIRS2 sounder uses a two field of view 
single cloud layer cloud filtering algorithm (i.e., equation 3, 4, and 5 in 
the text). Based on this study we found that the multiple layer cloud fil- 
tering algorithm showed significant improvement over the single cloud layer 
system. Therefore we are currently in the process of modifying our operational 
system to test whether this result holds true for real data analysis. 
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Figure 1. RMS temperature errors for AMTS and HIRS sounders (physical retrieval 
scheme and regression). 


74 



PRESSURE (mb) 


100 
114 
129 
147 
167 
190 
216 
245 
278 
316 
359 
408 
464 
627 
599 
681 
774 
880 
1000 
0.4 


Figure 2 , Compar 
errors 





PRESSURE (mb) 




RMS TEMPERATURE ERROR 


C»p.r1sor of retnovols “5 

cloS'«t"Jln“‘ra‘d5anc«'nJated at GLAS (RMS tenperatoro errors). 



7 



ORIGINAL PAC2 E9 
OF POOR QUALITY 


WATER VAPOR PROFILE RETRIEVALS FROM THE HIRS/MSU SOUNDER; CLEAR AND CLOUDY 

SIMULATION STUDIES 

D. Reuter and J. Susskind 


Currently, the GLAS physical retrieval scheme is capable of producing 
accurate atmospheric and surface temperature retrievals and accurate measures 
of cloud fields from radiances measured by the HIRS/MSU sounding unit (Susskind 
et al . . 1982). Since the presence of water vapor could contaminate these 
radiances only those infrared and microwave channels which are quite insensitive 
to the amount of water vapor are used to determine the temperature and cloud 
fields in the analysis. However water vapor is, in itself, a variable of 
meteorological interest and the HIRS instrument has a number of channels which 
are quite sensitive to the presence of water vapor; specifically the surface 
and lower tropospheric channels 8 and 10 and the mid and upper tropospheric 
channels 11 and 12. In order to determine if these channels contain any useful 
water vapor information we have carried out studies in which water vapor 
retrievals were performed, based on simulated clear and cloudy radiances and 
retrieved atmospheric and surface temperatures. The method used for simulating 
the clear and clody radiances has been described previously in this report 
(Reuter et al .) and thus no further detail will be given here. 

The water vapor retrieval scheme is based on considerations similar to 
those used in the retrieval of temperature information. Specifically, for two 
atmospheric profiles which have the same temperature structure, but whose water 
vapor contents are different by a small amount 6y(P) it may be shown that: 

T 

Tilii(P) +6 m(P)] - Tt[y(P)] =:/ Wi(P) 6 m (P) dtnp (1) 

Ps 

where T^ is the equivalent black body temperature for radiance R^ in channel i, 

T is the satellite pressure, Pg is the surface pressure and W'j is a water 
vapor sensitivity function for channel i given as 


Wi(P) 


oinp 

■■■■ ■ 


u(P). 


In equation 2 Ti-(P,e) is the frequency averaged transmittance in channel i 
from pressure P to 7 at angle e. 


( 2 ) 


In this very preliminary work two somewhat drastic, but not unreasonable 
assumptions were made to simplify the problem. The first assumption arises 
because each channels is primarily sensitive to changes to water in specific 
pressure intervals and thus what is determined j^n not 5 m (P) but rather a mean 

P 


layer change 6. The second assumption is that / Wj(P) dinp is essentially a 

— P$ 

constant W for each channel for the profiles studied. Under these assumptions 
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equation [1} becomes: 
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TiCu(P) + 6 m(P)] - TiCu(P)] = W 6. (3) 


This equatl0|n suggests that a regression treatment may be used to find ¥ as a 
function of I’ and brigtness temperature differences. The method of analysis 
Is, then, as follows: 

1) Simulated radiances are used as Input to the GLAS temperature 

scheme and retrieved temperature profiles are produced, 

2) These retrieved temperature profiles are used to simulate brightness 

temperatures In channels 8, 10, 11, and 12 using climatological 

water vapor profiles 

3) The data Is divided Into two sets, the first to generate the regression 
relationship and the second set to use as verification. The regression 
equation Is then solved to find the regression matrix A. 

(„TRUE . „CUM) . ;^(X^0BS . t^CALC) (4) 


where 

A = (aTT at + y)- 1 Ay ATT (5) 

and AT = Ay = - pCLIMj ^ is a ridge parameter. 

4) This regression relation is used on the second set of data to generate 
retrieval statistics. 

Note that the retrieved temperatures are used in the calculation of 
Equations (1) thru (3) are most nearly true if there is.no difference In 
the temperature profiles giving rise to T^^“^ and However, In the 

processing of real data, we would only have the retrieved temperature profiles, 
thus to attempt an accurate assessment of the ability of the HIRS Instrument to 
retrieve water vapor profiles, we have Included this source of noise In the 
simulation procedure. 

The results of this study are given in figure 1 in terms of the percent 
error of the retrieved water vapor column densities in the pressure regions 
1000-700 mb, 700-500 mb and 500-300 mb. Also given in this figure are the 
percent errors of the climatological guess. Results for two cases are shown: 
case $1 is for a very moist summer ocean tropical zone, and case Wi Is for 
a dry winter land mid-latitude zone. For case results are presented for 
simulated clear and cloudy conditions. Note that the presence of clouds does 
not degrade the results significantly. The lines are only used to connect the 
points and do not indicate accuracies In the regions between points. Table 1 
presents the absolute accuracies of the retrievals in terms of rms column 
densities for the pressure levels specified above. From this table It may be 
seen that even though the precent errors for Wi are much worse than those for 
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$ 1 , the absolute accuracies are somewhat better. 

Based on the results of this study it appears that there is usable water 
vapor information contained in the HIRS water vapor channels. We are presently 
in the process of refining this technique in the hope that we may impliment 
water vapor retrievals as part of the real data analysis scheme. 


Table 1 

RMS errors of retreived water vapor column densities (g/cm^) (independent set) 


P(mb) 

Si(clear) 

$ 1 (cloudy) 

Wi(clear) 

300-500 

0.029 

0.028 

0.012 

500-700 

0.14 

0.16 

0.063 

700-1000 

0.50 

0.58 

0.41 

1-1000 

0.51 

0.65 

0.43 
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Figure 1. Percent errors of water vapor retrievals and climatologies: points 
lie in the middle of the mean pressure regions. The lowest point 
is extended vertically to 1000 mb. 
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A LEVEL 2.4 SECOND-ORDER CLOSURE MODEL FOR THE PREDICTION OF TURBULENCE 

H. M. Helfand and J. C, Labraga 


Mellon and Yamada (1974) and Yamada (1977) have developed a hierarchy of 
second-order closure models for the prediction of atmospheric turbulence. 

These models range from the simple level 2 model, which consists of a system of 
purely algebraic eguatlons, to the more complicated, expensive and accurate 
level 4 model, which consists of ten prognostic equations for the second order 
turbulent moments. 

Yamada (1977) has proposed a level 2.5 model as a good compromise between 
the simplicity of the level 2 model and the accuracy of the level 4 model. 

This model has only one prognostic equation, that for the turbulent kinetic 
energy. 


1/2 q2 = 1/2 (u? + v^ + W^). 


( 1 ) 


We have analyzed the model, however, and we have found that It Is Ill-defined 
and singular when It predicts a value of 1/2 q2 significantly '.'mailer than that 
predicted by the level 2 model. Accordingly, we have created a new second- 
order closure model which combines the best features of the level 2 model and 
the level 2.5 model. We have labeled this the level 2.4 model. 

Both the level 2 and the level 2.5 models define turbulent transport coef- 
ficients for momentum and heat according to 


Kp, = qtS,n (2) 

and 

Kfj = qtS|^, (3) 

where q Is the turbulent velocity scale defined In (1), t Is a turbulent length 
scale, and S^, and represent the dimensionless parts of the coefficients. 

Sfli and Sn are functions of the dimensionless shear, 


S = t/q |Uz|. (4) 

and of the dimensionless Brunt-Va1sa11a frequency* 


M = t/q (g02/0o)l/2 , 

One can also express the anisotropy, 

V w^(l/3 q2) . (6) 

as a function of S and M. Is the ratio of the actual turbulent kinetic energy 


■••'■111 A 
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in the vertical mode to that energy which would be in the vertical mode if the 
flow were Isotropic. Physical reality demands that should be well within 
the limits, 

0 < Cy^ < 3. (7) 

Isotropy occurs when = 1. 

The level 2 model determines q as a function of the mean state of the atmo- 
sphere according to 

q2 = Bi *2 (1^ . Ri s,,) , (8) 


where Ri is the gradient Richardson number, 

R1 = • (’) 


and where Bj is a constant. Eq. (8) tells us that S and M are functions of Ri 

f>* r>j 

in the level 2 model. Sfy and therefore are also functions of Ri , and, 
as one can see in Fig. 1, they are well behaved. Cy^, for instance, is well 
within the bounds of (7) and it decreases slowly with increasing atmospheric 
stability. 

The level 2.5 model has an extra degree of freedom in comparison to the 
level 2 model. In the level 2.5 model, q is a prognostic variable and S and M 

mist now be expressed as functions of Ri and q. Thus S^, and Cy^, which 
depend on S and M exactly as they do in the level 2 model, are also functions 
of both Ri and q. 

While Eq. (8) does not necessarily hold for the level 2.5 model, it does 
hold whenever q reaches its equilibrium value of qe(Ri)» Thus, the level 2 model 
is a special case of the level 2.5 model, and Fig. 1 represents the Ri dependence 

of Sm, Sf, and Cy^ for the level 2.5 model for the special case q = qe (Ri)« 

Fig. 2 represents the behavior of these functions for the case q = 1.4 qe (Ri) 
or the case where the turbulent kinetic energy is double its equilibrium value. 

Here Sm, Sj, and Cy^ behave quite smoothly, and, in fact, are less sensitive to 
Ri than they were in the level 2 model. 

When q = .7 qe (Ri) or the turbulent kinetic energy is half of its equilibrium 

^ nj 

value, on the other hand, the behavior of S^, Sj, and Cy^ is singular. Fig. 3 
demonstrates that these functions approach plus infinity as the Richardson number 
decreases to -.68 and then they increase from minus infinity to negative finite 
values. This behavior is nonphysical, and the condition (7) is clearly violated 
for Ri < -.55. The level 2.5 model is basically an equilibrium model for turbulence 
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but When 1/2 q2 is significantly smaller than Its equilibrium value one cannot 
assume that the other turbulent quantities are in an equilibrium state. 


To remedy the problem with the level 2.5 model for q < qe» we have proposed 
the following model which we term the level 2.4 model. Define and 
according to (2) and (3) but take 


^m “ ^m ) 
Sh • Sh («i) 

from the level 2 model for q < qe (Ri)» 4 nd take 


( 10 ) 


^m * ^m (*^^ » ^) 
Sh = Sh (Ri, q) 


( 11 ) 


from the level 2.5 model for q >. ^le (Ri)* The value of q in (4) and (5) is the 
prognostic one in either case. This model has the advantage of allowing and 
Kh to be proportional to a prognostically determined turbulent velocity scale 
yet it assures one that the Richardson number dependence of these coefficients 
will be reasonable. It also allows turbulence to exist when Ri exceeds its 
critical value, as does the level 2.5 model, for in that case, q^ (Ri) = 0 and 
q > 0 Implies that the model is identical to the level 2.5 model. 
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model when the turbulent kinetic energy is double Its equilibrium value. 
Notation is that of Fig. 1. 
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Dependence of S^, S>, and on Richardson number for the level 2.5 
model when the turbulent kinetic enerqy is half its equilibrium value 
Notation is that of Fiq. 1. 


SELECTION OF A BEST CANDIDATE HIGHER-ORDER CLOSURE SCHEME 
FOR TURBULENCE IN THE GLAS FOURTH ORDER GCM 

J. C. Labraga and H. M. Helfand 


Mel lor and Yamada (1974) have presented a hierarchy of second-order closure 
models for turbulence ranging from the complete level 4 model, with ten prognostic 
equations, to the simple level 2 model, which is a fully algebraic system of 
equations. On the basis of a set of Wangara simulations, they concluded that 
an intermediate member of the hierarchy, the level 3 model with two prognostic 
equations and eight diagnostic equations, was the optimum compromise between 
expense and accuracy for geophysical applications. 

Later, Yamada (1977) introduced a yet simpler compromise between the level 
3 and the level 2 models. His level 2.5 model carried turbulent kinetic energy 
as its only prognostic variable with the remaining seven variables determined 
prognostical ly. Miyakoda and Sirutis (1977) have selected this level 2.5 model 
to represent the vertical turbulent fluxes of heat, momentum and moisture in 
the planetary boundary layer (PBL) and in the free atmosphere in the GFDL 
General Circulation Model (GCM). 

Elsewhere in this volume, Helfand and Labraga have demonstrated that the 
level 2.5 model is ill-defined and singular. In its stead, they have developed 
a level 2.4 model which incorporates the best features of the level 2.5 and 
level 2 models and yet is not singular. 

We have tested the level 2.4 model by repeating a 48 hour simulation of 
the Wangara experiment carried out by Yamada and Mellor (1975) with the level 3 
model. Our model had 31 layers progressively concentrated near the ground so 
as to have vertical resolution similar to that used by Yamada and Mellor. We 
also repeated our simulation with a simpler level 2 model so as to see where 
the level 2.4 simulation fell between those by its simpler and mort- complicated 
cousins. 

We discovered that the level 2.4 model gave results remarkably similar to 
those of the level 3 model (see, for example. Fig. 1). We also saw, to our 
surprise, that the simulation by the level 2 model was virtually identical to 
that by the level 2.4 model. 

To examine the feasibility of using the level 2.4 scheme in the GLAS Fourth 
Order GCM (Kalnay-Ri vas et al . , 1977), we repeated the Wangara simulation at a 
drastically reduced vertical resolution. In this experiment, the PBL was 
resolved by only four grid points, only three of which were located within 2 km 
of the ground. The layers were chosen as those of the PBL-simulating version 
of the GLAS Fourth Order Model. The surface flux parameterization was also 
taken from that model. This parameterization uses an extended Monin-Obukhov 
similarity scheme to predict surface fluxes on the basis of the temperature and 
wind gradients over the lowest 150 m of the atmosphere. 

Fig. 2 shows that even at this drastically reduced vertical resolution, 
the level 2.4 model gives a good qualitative simulation of the diurnal cycle. 

Fig. 3 shows that the level 2.4 model, together with the surface flux parameteri- 
zation does an excellent job at predicting surface fluxes at the coarse resolution. 
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The results of the experiment demonstrate the suitability of the level 
2,4 model together with the surface flux scheme for use in the GLAS model. 

The experiment was repeated with the level 2 model, and the results were 
again indistinguishable from those of the level 2.4 model. The level 2 model is 
cheaper to run than the level 2.4 model because it requires fewer computations, 
it is less subject to numerical instability, and it does not require storage of 
a turbulent kinetic energy field for prognostic purposes. Its accuracy is 
apparently identical to that of the level 2.4 model, and therefore, on the 
basis of the present study, the level 2 model is the best candidate for use in 
the GLAS Fourth Order GCM. 
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FIGURE CAPTIONS 


Fig. 1, Time and space variation of the virtual potential temperature 0y in 
units of degrees C for the observations (a) and for the simulation 
by the level 3 model (b) and in units of degrees K for the simulation 
by the level 2.4 model (c) and by the level 2 model (d). 


Fig. 2. Virtual potential temperature Gw profiles at 1200, day 33 to 0900, 

day 34 from the observations (a), from the simulation by the 31-layer 
model (b) and from the simulation by the 4-layer model (c). 


Fig. 3. Hourly variation of the surface heat flux from observations (triangles), 
from the simulation by the 31-layer model (solid line) and from the 
simulation by the 4-layer model (dashed line). The units are m s"l K. 
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A STATISTICAL DYNAMICAL STUDY OF LARGE-SCALE FLOW ANOMALIES 

S. 0. Schubert 


1. INTRODUCTION 


In this study the large-scale intraseasonal variability of the winter time 
extratropical atmosphere is examined within the framework of a hemispheric 
equivalent barotropic model expanded in horizontal empirical orthogonal functions 
(EOFs) in the wave domain. This approach allows for a detailed observational 
and statistical analysis of the anomaly fields and in addition provides insights 
into the dynamical processes governing their behavior. 


2, OBSERVATIONS 


The anomaly fields are computed as deviations from a mean seasonal cycle 
based on a ten-year data set of the 500 mb streamfunction. A statistical test 
developed by Preisendorfer et al» (1981) shows that the first 12 time dependent 
principal components (PCs) associated with the EOF anomalies have a variability 
which is significantly different from that of PCs computed from spatially and 
temporally uncorrelated noise. The spatial patterns of the statistically signi- 
ficant modes of variability are examined in the context of the time mean kinetic 
energy spectrum and are interpreted by superimposing them on the winter mean 
flow. 


The first EOF is dominated by zonal wave numbers 0 and 2 and has a spatial 
pattern resembling the Pacific/North American pattern reviewed recently by 
Wallace and Gutzler, 1981. When superimposed on the mean flow with large positive 
and negative PCs (2.5 standard deviations) this mode exhibits a behavior 
resembling an index cycle in the western hemisphere. 

The second EOF has a large anomaly center over the North Atlantic and is 
dominated by zonal wave numbers 1 and 2. This mode is associated with fluctua- 
tions in the North American east coast jet and seems to be related to the North 
Atlantic Oscillation also reviewed recently by Wallace and Gutzler, 1981. 

The third EOF pattern is primarily composed of zonal wave numbers 1 and 3. 
This mode has anomaly centers closely coupled to the positions of the east 
coast jets and when superimposed on the mean flow shows, for the positive case, 
a very zonal flow over the U.S. and for the negative case a high amplitude wave 
pattern with a strong North American east coast jet and extensive Pacific 
trough. This EOF also bears some resemblance to the Pacific/North American 
pattern. 

A comparison with the occurrences of blocking in the Pacific and Atlantic 
as determined by Hartmann and Ghan, 1980 suggests that the fourth and ninth 
EOFs are related to Pacific and Atlantic blocking, respectively. The fourth 
EOF is primarily associated with zonal wave numbers 1 and 3 and has a large 
anomaly in the polar regions extending southward over the west coast of North 
America, Greenland and Asia. The spectral defcomposition and spatial pattern 
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for EOF nine 1s shown In Flqs. la and lb, respectively. This nwjde its § large 
anomaly In the eastern North AtlantlG with an afloi?ialy os opposite sign to the 
south. It Is dominated by zonal wave numbers 2 and 3 ind a meridional scale 
corresponding to approxltnately 1,5 wavelengths between the equator and North 
Pole. Here n-m Is Is the number of zeros the spherical hartiwnlc has between 
the poles. 


3. MODELING 


The equivalent barotropic model Is formulated as a prognostic equation for 
the anomaly streamfunctlon PCs at the 500 mb level. The unknown constants 
Involved In the boundary layer parameterizatlons are determined by regression 
methods. The seasonal cycle enters Into the anomaly equations via the advec- 
tlve terms and a san additional source term due to the fact that the mean terms 
do not by themselves satisfy the vorticity equation. 

The barotropic stability of the zonal 1y varying winter mean flow may be 
tested with respect to perturbations composed of the EOF patterns. For this 
analysis, one PC is allowed Into the model at a time to determine whether such 
a perturbation will grow or decay by gaining or losing energy from the mean 
state. 


Fiq. 2 shows the growth rates for the first 25 PCs without a dissipative 
term In the model. These results suggest that, except for the 8th PC, all the 
statistically significant modes have the potential to grow via kinetic energy 
transfers from the mean flow. The second and third PCs are the most unstable 
having e-folding times of a little more than two weeks. An analysis of the 500 mb 
winter mean flow shows that the zonally averaged state Is barotropically stable. 
Thus the zonally asymmetric components of the mean state which account for less 
than 16 % of the mean flow kinetic energy and are composed primarily of zonal 
wave numbers 1,2 and 3 are key factors influencing the stability of the EOF 
patterns. 

A more complete normal mode analysis shows that with 40 PCs the fastest 
growing mode is also fixed in space and is composed primarily of EOFs 4 and 12. 
However, it should be noted that there is in general a strong sensitivity of 
the growth rates and period of the normal modes to the level of EOF truncation. 
These results are qualitatively very similar to the results obtained for the 
climatological mean January 300 mb flow by Simmons et a1 . , 1982, In their study 
of barotropic wave propagation and instability. They found the fastest growing 
mode to have a period of about 50 days and an e-folding time of about one week. 

By varying the base state they also obtained some modes with Infinite period. 

The nonlinear model is analyzed in the context of steady state solutions 
employing a high degree of EOF truncation. The EOF interaction coefficients 
are linear combinations of the usual spectral Interaction coefficients. An 
examination of these coefficients shows that by far, the strongest interactions 
Involve the 9th EOF. 

Multiple equilibria are found for severe truncations (less than eight modes) 
without including orography In the model. In this case, the advectlon of the 
winter mean flow vorticity by the anomalies plays the role of an orographic type 
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term. This term was found to dominate the topographic effects at the 500 mb 
level. 


4. CONCLUSIONS 


A statistical and observational analysis of the EOFs computed from a ten- 
year winter data set has shown that as many as the first 12 patterns may be 
associated with physically significant modes of variability. The first three 
EOFs show some relationships to index cycles and/or some of the teleconnection 
patterns discussed recently by Wallace and Gutzleri 1981. The fourth EOF Is 
related to North Pacific blocking and the ninth bears a striking resemblance to 
the classical North Atlantic block. Details of the EOF patterns must, however, 
he viewed with some caution since the limited extent of the data set allows 
only a marginal resolution of the "true" climatological variances and covariances. 

The linear normal mode analysis of the EOF model shows that all of the 
statistically significant EOFs except for the 8th have the potential to grow at 
the expense of the mean flow. This supports the recent work by Simmons et al. . 
1982, which also suggested that a major source of the observed variability of 
the large scale extratropical winter time circulation Is derived through baro- 
troplc Instability of the zonal ly varying mean flow. 

Multiple equilibria have been obtained for highly truncated versions of 
the EOF model In which orographic effects are neglected. In these cases, the 
advectlon of mean flow vortlclty plays the role of an orographic type term. An 
examination of the EOF Interactions shows that the dominant Interactions Involve 
the North Atlantic blocking mode (EOF 9). 
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NORMAL MODES OF THE GLAS 4th ORDER MODEL 
S. Bloom 

ORIGINAL PAGE IS' 

1. INTRODUCTION OP POOR QUALITY 


A recurring problctn in numerical weather prediction is the initialization 
of primitive equation models. Recent advances in this area have involved the 
development of nonlinear normal mode initialization techniques (Machenautr, 
1977; Temperton and Williamson, 1981; Williamson and Temperton, 1981). The 
computation of the normal modes of the (linearized) GLAS 4th order model is the 
first step in an investigation of the impact of these initialization techniques 
upon the full GLAS model. This note discusses some of the features of the GLAS 
model's normal modes. 


2. GLAS MODEL MODES 


a) Vertical Structure 

The ^LAS_4th order model equations are linearized about a basic state 
at rest (u = v = 0) and with an assumed average vertical temperature profile, 
T(o), taken from the ICAO standard atmosphere. Using methods very similar 
to those described by Temperton and Williamson (1981), the vertical structure 
of the model modes is separated out. This procedure involves combining the 
linearized thermodynamic and hydrostatic equations in their finite differenced 
form. Because the vertical and horizontal structure is assumed to be uncoupled, 
a vertical structure problem is obtained. The solution to this vertical struc- 
ture equation is an eigenvalue problem. Since the version of the model used 
here has nine o levels, there are nine eigenvalues/eigenvectors, each corre- 
sponding to a particular vertical mode. When the linearized GLAS model equa- 
tions are projected onto these vertical modes, they are transformed into the 
linear shallow water equations in spherical coordinates: 


3U. 

3t 

- 

.a...-., yiA. = 0 

a cos ax 

(la) 


+ fu, + 

-i- “ 

(lb) 

3hn 

at 

a cos <|i 

^ (V,cos ♦)] = 0 

(Ic) 


Here Dj^ is the tth eigenvalue and (U)i, vji, he) are the variables 
associated with the tth vertical mode. Note, On is often referred to as an 
"equivalent depth". 
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b) Horizontal Structure ’ 

Following Temperton (1977) and Dickinson and Williamson (1972 )» the fields 
are expanded in zonal Fourier harmonics, and scaled: 



so that eq. (1) can be written for each latitude (j>j: 


(-i ) cos (|»j au* 

; at j 

(-i ) cos (j>j • av* 

V at , j 


- fj cos = 0 

- f . cos d> .uij - cos ^ . f . = 0 

J J J a V 7 J 


(2a) 

(2b) 


(-1) cos (j»j ah*^ 

V at /j 


+ s *Cg y* + c# a 

a J 


(v*cos a). = 0 

a a(j> J 


(2c) 


with Cji = (g0j^)l/2 and S' = sin(SAX) '4 _ cos(SM)l‘ 

3ax ' i-l 


result of 4th order differencing in x. If the modes have a harmonic time depen- 
dence (i.e., e^''^), then for a sufficiently small At: (-i) a / y. = » ^( )• 

With the use of 4th order differencing in <)>, eq. (2) can be expressed as a matrix 
eigenvalue problem: 


L u - V u = 0 


( 3 ) 


where 


“j = 



* (cos 


and the matrix of coefficients L is real and symmetric. Thus for a given zonal 

M 

wavenumber s, the eignevectors of L give the latitudinal structure of the hori- 

M 

zontal modes (which have eigenfrequencies v). 
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c) Degrees of Freedom - Influence of Poles OF POOR Ol’AHTY 

The number N of eigenvectors and eigenvalues in eq. (3) reflects the number 
of degrees of freedom in the GLAS model for a given zonal wavenumber. Since 
the horizontal modes can be divided into two classes (symmetric or antisymmetric 
across the equator), they need to be computed only in one hemisphere. 23 latitude 
points in a model with a meridional spacing of 4®. Eq. (2) suggests that N = 

69. but consistency requirements at the pole reduce N: 


S = 0 ; u* = V* =0 and h* 0 (N = 67) 

pole pole pole 

S = 1 ; h* , =0 and u* = + v* , 0 (N » 68) 

pole pole “pole 

{+ North. “ South) 

S > 2 ; u* = V* = h* = 0 (N <= 66) 

“ pole pole pole 

A problem arises with the S = 1 modes, due to the use of centered differencing 
on a nonstaggered grid in the GLAS model. This can be seen by comparing the 

roles of h* and v* in eq. (2); where eq. (2b) relates v* to h* ^ and h* . but 

pole 86® 82° 

V* does not influence h* „ or h* „ in eq, (2c). Thus the matrix L is not 
pole 86® 82® ^ ' 

symmetric for S = 1. The degrees of freedom at the pole are represented by 
inertial oscillations (one physical, one computational) which are completely 
decoupled from the remaining modes. 


3. COMPARISON OF GLAS MODES WITH HOUGH MODES 


Longuet-Higgins (1968) and Kashara (1976) discuss alternative ways to 
extract the horizontal structure from eq. (1) through the use of spherical 
harmonic expansions. The resulting horizontal modes are called Hough modes, and 
they represent an analytic solution to eq. (1) to which we can compare the GLAS 
modei modes (and eigenfrequencies). Using an algorithm similar to the one used 
by Kasahara (1976), Hough modes with the same Dji's as the GLAS modes were 
computed. Both the GLAS and the Hough modes can be labeled by the number 
of their latitudinal nodes; i.e. a "low" mode has little structur<i, while a 
"high" mode has many oscillations. What follows are some preliminary results 
of the comparisons between GLAS and Hough modes. 

(i) The GLAS modes consist of two types: physical and computational. The 
computational modes are an artifact of the finite differencing approximations 
and they do not represent true solutions to eq. (1). There should be roughly 
the same number of computational and physical modes, although identifying the 
higher physical GLAS modes becomes increasingly more difficult as their struc- 
ture increases. 
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(ii) For low wavenumbers and the lowest modes, there is eitcellent agree- 
ment between the etgenfrequencies of the Hough and 6LAS modes. The agreement 
deteriorates as one compares the eigenfrequencies of progressively higher modes. 

(iii) An examination of the polar and near-polar structure of the lowest 
modes for S = 0, 1 and 2 shows the nature of the GLAS model's (linear) treat- 
ment of the pole. The best correspondance between Hough modes and GLAS modes 
near the pole occurs for S = 2, because no polar approximation is required for 
the computation of the GLAS modes for zonal wavenumbers greater than one. 

The greatest discrepancies near the pole occur for S = 1; this is consistent 
With the model's difficulty in handling the polar wind. There is much better • 
agreement between GLAS and Hough modes near the equator for S = 1. Finally 
for S = 0, the agreement near the pole lies in between those for S = 2 and 
S = 1. This occurs because an approximate equation for h* is obtained by 

pole 

performing polar-cap averages of eq, (2c); this approximation somewhat degrades 
the GLAS modes near the pole for S = 0. 


4. SUMMARY 


A set of normal modes for the 4th order GLAS model has been computed and 
compared with a corresponding set of Hough modes. The agreement between the 
lower GLAS (physical) and Hough modes is a useful check on the accuracy of the 
GLAS mode computation. Future work with the GLAS modes will involve: the 
diagnosis of data, linear initialization experiments, and the development of a 
nonlinear normal mode initialization process for the full GLAS 4th order model. 
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LAGGED AVERAGE FORECASTING, AN ALTERNATIVE TO MONTE CARLO FORECASTING 

R. N. Hoffman and E, kalnay 


In order to use the information present In past observations and simul- 
taneously to take advantage of the benefits of stochastic dynamic prediction 
we formulate the lagged average forecast (LAF) method. Jn a LAF, just as in a 
Monte Carlo forecast (MCF), sample statistics are calculated from an ensemble 
of forecasts. Each LAF ensemble member is an ordinary dynamical forecast (ODF) 
started from the initial conditions observed at a time lagging the start of 
the forecast period by a different amount. These forecasts are averaged at 
their proper verification times to obtain a LAF. The LAF method is opera- 
tionally feasible since the LAF ensemble members are produced during the normal 
operational cycle. 

To test the LAF method, we use a two layer, f-plane, highly truncated 
spectral model, forced by asymmetric Newtonian heating of the lower layer. In 
the experiments a long run is generated by the primitive equation version of 
the model which is taken to represent nature, while forecasts are made by the 
quasi geostrophic version of the model. On the basis of forecast skill the LAF 
and MCF are superior to the ODF; this occurs principally because ensemble 
averaging hedges the LAF and MCF toward the climate mean. The LAF, MCF and ODF 
are all improved when tempered by a simple regression filter; this procedure 
yields different weights for the different members of the LAF ensemble. The 
tempered LAF is the most skillful of the forecast methods tested. The LAF and 
MCF can provide a priori estimates of forecast skill because there is a strong 
correlation between the dispersion of the ensemble and the loss of predictabi- 
lity. In this way the time at which individual for.tjcasts lose their skill 
can be predicted (Fig. 1^. This is not easy since skill of individual fore- 
casts, as measured by D, ',n RMS measure of forciu?st "kill, breaks down suddenly 
anywhere from four to twenty days into the forecast (Fig. 2). 

These results and the application of the LAF method to more realistic 
models and to monthly or seasonally averaged forecasts are discussed elsewhere 
in more detail (Hoffman and Kalnay, 1983). 
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Fig. 1. Scatter plots of predicted versus observed time for D to reach a value of 
0.5 for the LAF (a) and the MCF (b). 1) is the ensemble rms value of D. 
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Fig. 2. 


Time evolution of D for the ODF experiments. Shown 
are selected portions of the evolution of D for each 

case (light ■■ - ■). These curves have been smoothed 

to eliminate the roughness caused by white noise 

observing error. Also shown are <D^> (•— ) 

>)• Error bars are calculated 
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A STOCHASTIC-DYNAMIC MODEL FOR THE SPATIAL STRUCTURE OF FORECAST ERROR STATISTICS 
R. Balgovind, A. Dalcher, M, Ghtl and E. Kalnay 


A simple model that yields the spatial correlation structure of global 
atmospheric mass-field forecast errors is derived. The model states that the 
relative potential vorticity of the forecast error is forced by spatially 
multi -dimensional white noise. The forecast error equation contains a non- 
dimensional parameter Cq which depends on the Rossby radius of deformation. 

From this stochastic-dynamic equation, a deterministic equation for the spatial 
covariance function of the 500 mb geupotential error field is obtained. 

Three methods of solution are examined: 1) an analytic method based on 

spherical harmonics, 2) a numerical method based on stratified sampling of 
Monte-Carlo realizations of the stochastic-dynamic equation, and 3) a combined 
analytic-numerical method based on two successive applications of a fast Possion 
solver to the deterministic covariance equation. The three methods are compared 
for accuracy and efficiency, and the third, the combined method is found to be 
clearly superior. 

The model's covariance function is compared with global correlation data 
of forecast-minus-observed geopotential fields for the DST-6 period February- 
March 1976. The data are based on the GLAS forecast-assimilation system In use 
at that time (Ghil et al ., 1979). 

The model correlations agree well with the latitude dependence of the data 
correlations. The fit between model and data confirms that the forecast error 
between 24h and 36h is largely random, rather than systematic; the value of the 
parameter Cq which gives the best fit suggests that much of this error can be 
attributed to baroclinic, rather than barotropic effects. Deterministic influ- 
ences not included in the model appear at lEh and 48h. They suggest possibili- 
ties of improving the forecast system by a better objective analysis and 
initialization procedure, and a better treatment of planetary-wave propagation, 
respectively. 

An analytic formula is obtained which approximates well locally the model's 
global correlations. This formula is convenient to use in the calculation of 
weighting coefficients for analysis and assimilation schemes. It shows that 
Gaussian functions are a poor approximation for the forecast error correlations 
of the mass field. The first and second derivatives of the Gaussian functions, 
conmonly used to model geostrophic wind error correlations, are shown to be 
even less accurate. 
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Figure 1, Comparison of model correlations with the data correlations obtained 
when using 48h GLAS forecasts and the corresponding NMC analyses. 
Model results have c° = 190; for this value (f Cq, the results 
were almost indistinguishable from those with Cq = 230. 
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Figure 2a. Analytic approx imati on of global model correlations. Fig. 2a shows 
the shape of the Bessel -function correlations derived in this paper 
for different values of a = o(e). 
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Figure 3. Comparison between analytic formulae used to approximate forecast 
error correlations. The four functions compared are the Gaussian, 
currently used In 01 (marked Gauss In the figure), the Bessel func- 
tion of Figure 2 (marked j^i), the exponential (exp), and the one- 
dimensional formula of tha analytic model ( 1»D ). Figure 3a shows 
the functions themselves. Figure 3b their f 1 rst derivatives, and 
Figure 3c the second derivatives. According to the geostrophic 
assumtplon of operational 01, if either curve In Figure 3a were 
used for (<|>-1») correlations. Figure 3b would give the corresponding 
(♦-v) correlations, and Figure 3c the (v-v) correlations, up to 
a multiplicative factor. Notice large differences between the 
curves In Figure 3b, c. 
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Figure 3b 
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HIGH LATITUDE FILTEKING IN GLOBAL GBID POINT NOOELS 
L. L. Takacs and R. C. Balgovind 


I. INTRODIIQTION 


Thn process of filtering in qlobal prid point nwdels Is fairly common in 
meteoroloqical studiosi High latitude filtering as described in Gates et al . 
(1071) is done in order to avoid the use of prohibitively short time steps 
required duo to fast moving inertia-gravity waves near the poles. By applying 
filtering to certain quantities in the governing equations, the use of larger 
time-steps to run the model may be obtained. The specific filtering response 
needed to ensure stability is defined by the space and time difference schemes 
used,, and by the type of grid structure in the model (Arakawa and Lamb, 1977), 

We have found that when certain techniques are used in the application of 
the hiqh-latitude filter, serious distortions of the solution occur at all 
latitudes, associated with spurious energy and momentum transfers. These 
errors are a result of the inability of the specific filtering technique used 
to maintain in finite-difference form the irrotationality of the gradients of 
the streamfunction and velocity potential. 

In this paper we re-examine the effects of the distortions caused by 
filtering, and study its relationahip to various filtering techniques. In 
Section 2, a description of the model used for this study is presented, while 
Section 3 compares iiK)deT results. Conclusions are presented in Section 4 . 


2. MODEL DESCIUPTION 


The model used to examine the effects of filtering is based on the shallow 
water equations on a sphere. All variables are defined on a non-staggered 'A' 
grid (Arakawa, 1972) having a resolution of 10 degrees in longitude and 8.18 
degrees in latitude. Spatial derivatives are approxinmted by energeticaTly- 
conservative fourth-order differences (Kalna.y-Ri vas et al. . 1977), while the 
time derivatives are approximated by the leapfrog scliemeT” 

As the initial condition for these experiments, we followed Phillips 
(1959) by taking an initially non-divergent velocity field together with a 
height field that satisfied the balance equation, (see Pig. 1). This initial 
condition produces a super-rotation of the atmosphere In addition to a Rossby- 
Haurwitz wave which moves westward relative to the fluid. For a non-divergent 
barotropic atmosphere, the total solution will move from west to eadt without 
change of shape. Since we are using a free upper surface, however, some 
changes in the solution can be expected due to the divergent nature of our 
model , 

For our experiments, we compare the effects of four filtering techniques. 
The first technique, as employed by Kalnay-Rivas et al . (1977), uses Fourier 
filtering of thp prognostic variables u, v and h at every time-step. This 
technique does not conserve total energy nor does it preserve the irrotatlonal 
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properties of the system. Since this fIVter Is applied to the prognostic fields 
themselves, the effects of the filter are cumulative In time. This technique 
will be referred to as technique (b) In future reference. 

The second technique used for filtering Is similar to the first except 
that the time-changes of the prognostic flux and mass fields rather than the 
fields themselves are filtered, I.e., h’u, hv, and fi, where (*) refers to the 
time-change of the quantity. Smoothing of the time derivatives to prevent 
linear computalonal Instability has been employed by Vanderman (1972) and 
Temperton (1977). This technique also suffers from non-conservation of the 
energetics and Irrotatlonnllty of the system, although the effects of the 
filter are non-cumulative since only time-changes are smoothed. This technique 
will be referred to as technique (c) In future reference. 

The third technique follows that of Arakawa and Lamb (1977). In this 
technique, the zonal pressure gradient and zonal mass flux terms are filtered 
In such a way that the globally Integrated total energy of the system Is con- 
served. Although this technique was originally designed for an energy and 
enstrophy conserving scheme on a 'C* grid, our Interest here Is to study Its 
performance In relation to our fourth-order energy conserving scheme on the ‘A* 
grid. This technique will be referred to as technique (d) In future reference. 

Thf; last technique, which will be denoted as the P.6, (for pressure gra- 
dient) filter. Is one in which the vorticity of the system Is conserved In 
finite difference form. In a personal communication. Dr. Akio Arakawa pointed 
out that filtering the prognostic variables, their tendencies or the zonal 
pressure gradient and mass flux terms may not preserve the Irrotatlonallty of 
the pressure gradient In finite difference form due to the latitudinal depen- 
dence of the filter response. Because of this, a computational source of vor- 
ticity Is generated near the poles which ultimately affects the solution at 
all latitudes. By a priori filtering of the heights used In the pressure 
gradient terms In tWe momentum equations, the Irrotatlonal properties of the 
system are preserved and no spurious energy or momentum transfers results. 

This technique will be referred to as technique (e) In future reference. 


3. RESULTS 


Figs. 2, 3 and 4 show the time evolution of the height field at days 4, 8 
and 12. In each figure, five plots are shown corresponding to the four tech- 
niques discussed in Section 2 in addition to a non-filtered control run. The 
non-flltered control run Is taken to be truth In these experiments, and will 
be referred to as technique (a) in future reference. 

We can see from these figures that all techniques perform reasonable well 
through the first four days of Integration. By day 8, however, techniques (b) 
and (d) show major dissimilarities compared to techniques (a), (c) and (e). It 
should be noted that the control run has mo ’^d approximately 90 degrees eastward 
as expected, with very little change In shape. In techniques (b) and (d), the 
amplitude of the Rossby-Haurwitz wave has decreased considerably. Also, technique 
(c) shows tilting east of the trough which Is not apparent In the control run nor 
In technique (e). Finally, by day 12, only technique (e), which conserves vor- 
ticity In finite difference form, shows success In' simulating the control run. 
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Fiq. 5 shows the time evolution of the globally Integrated kinetic energy 
due to the eddy flow. We see from the control run that the eddy kinetic energy 
decreases slightly from Its Initial value during the 12-day run. There Is also 
a corresponding Increase In the kinetic energy due to the zonal flow (not shown) 
which maintains total energy conservation. Techniques (b), (c) and (d) show 
significant deviations from the control run which are a result of the non-con- 
servation of vorticity Inherent In those techniques. Physically, the change 
In vorticity near the poles generates tilting of the trough-ridge lines which 
result In strong equatorward transports of u-momentum at mid-latitudes. This 
equatorward transport triggers the barotropic Instability mechanism which pro- 
duces a transfer of eddy kinetic energy to zonal kinetic energy. However, It 
can be seen that the P.G. filter, technique (e), produces almost Identical 
results compared with the non-flltered control run. 


4. CONCLUDING REMARKS 


We have seen that certain high-latitude filtering techniques create 
spurious energy and momentum transfers, affecting the solution at all latitudes. 
These spurious transfers are due to a non-conservation of vorticity under the 
pressure gradient force, as well as to a non-vanishing curl In finite-difference 
form of the gradients of the streamfunctlon and velocity potential. By a priori 
filtering of the heights used In the pressure gradient terms In the momentum 
equations, the Irrotatlonal properties of the system are preserved and no 
spurious energy or momentum transfers result. This technique has shown super! ar 
results compared with filtering prognostic field variables, their tendencies or 
filtering with energy-conservative techniques* 
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LONGITUDE 


LONGITUDE 


Contoured height fields at day 4 using a) no filter, b) prognostic 
field filter, c) tendency filter, d) energy-conserving filter, e) P.6, 
filter. 
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Fiq. 5. Globally averaqed and normalized eddy kinetic energy for the a) no 

filter, b) prognostic field filter, c) tendency filter, d) energy-con- 
servinq filter, e) P.G. filter. 
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LARGE AMPLITUDE STATIONARY ROSSBY WAVES IN THE SOUTHERN HEMISPHERE: 

OBSERVATIONS AND THEORY 

E. Kalnay and J. PaegTe 


Studies by Van Loon and Jerine (1972), Van Loon et a1 . (1973), Trenberth 
(1979) and others indicate that stationary waves in the southern hemisphere are 
of planetary scale, and much weaker than those in the northern hemisphere, 

Kalnay and Halem (1981) reported the presence of large amplitude, short stationary 
waves during the month of January 1979 in the lee of South America and their 
disappearance in February 1979, In this paper we present further observational 
evidence of the January waves and discuss their possible origin. The principal 
results are: 

a) Large amplitude stationary Rossby waves with zonal wavenumber'^ 7 were 
present between 20‘’S and 40"5 both in the South Pacific and east of South 
America during January 1979, They appear in satellite observations as enhanced 
bands of high clouds associated with the South Pacific convergence zone and the 
Amazon, Examination of satellite observations during 1974-1979 indicates some 
correlation between the intensity of stationary cloud bands in the two regions, 

b) The stationary waves in the lee of South America are not of orographic 
origin since they are associated with a ridge rather than an trough east of the 
Andes, 

c) The waves could not be produced by the CISK mechanism suggested by Kalnay 
and Halem (1981), because of their essentially barotropic vertical structure, 

SST anomalies in the South Atlantic are of the same short scale as the waves 
during this period. However, the fact that they are stronger at the end of 
January^ and their phase relationship, indicated that they are a result and not 
a cause of the atmospheric waves. 

d) Tropical heating in the Pacific and South America as diagnosed by the 
6LAS FGGE analysis was stronger in January than in February 1979. At the same 
time the easterlies were much weaker, becoming westerlies in most of the western 
hemisphere. The January waves appear to have been produced by the enhanced 
tropical heating which produced local convergence of westerly momentum and 
allowed (Figs. 1 and 2) energy propagation to the subtropics unimpeded by the 
presence of a critical line. 

e) Numerical experiments with a 2-level hemispheric quasl-geostrophic 
model are in qualitative agreement with these conclusions. 
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A MODEL TO DETERMINE OPEN OR CLOSED CELLULAR CONVECTION 


H. M. Helfand and E. Kalnay 


A simple mechanism Is proposed to help explain the observed presence In 
the atmosphere of open or closed cellular convection. If convection Is produced 
by cooling concentrated near the top of the cloud layer, as In radiative cooling 
of stratus clouds, It develops strong descending currents which are compensated 
by weak ascent over most of the horizontal area, and closed cells result. 
Conversely, heating concentrated near the bottom of a layer, as when an air 
mass Is heated by warm water, results In strong ascending currents compensated 
by weak descent over most of the area, or open cells. This mechanism, unlike 
that of Hubert (1966), does not Invoke vertical variations of the eddy diffusion 
coefficients, and Is similar to the one suggested by Stommel (1962) to explain 
the smallness of the oceans' sinking regions. 

This mechanism Is studied numerically by means of a two-dimensional, non- 
linear Boussinesq model. For this purpose we define open (closed) convection 
as a convective circulation pattern In which the majority of the fluid has a 
descending (ascending) motion. An Internal heat source-sink destabllzes a 
layer of fluid adding no net heating. A steady state Is attained. The resulting 
circulation Is closed, as expected, when the cooling Is concentrated near the 
upper surface, and the heating Is spread throughout the Tower region. The mean 
lapse rate Is unstable In the upper half of the fluid and stable In Its lower 
half. Conversely, the circulation Is open when heating Is concentrated near 
the bottom. In this case, the lower half of the fluid has an unstable mean 
lapse rate, and the upper half of the fluid Is stable. 

The numerical results Indicate that the width of the plume produced by the 
cooling In the upper part of the layer or by the heating In the lower part of 
the layer Is largely Independent of the degree of vertical asymmetry of the 
Internal heating profile. On the other hand, compensating motion occupies 
a region which becomes broadi?«* as the heating profile becomes more asymmetric. 

In other words. If cooling Is very concentrated near the top of the layer with 
heating spread throughout the rest of the region or if heating is very concen- 
trated closed or open cells have an aspect ratio much larger than 1. These 
results may help explain the large aspect ratios observed in atmospheric con- 
vection. 
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UNSTABLE STEADY STATES IN QUASI-6E0STR0PHIC ATMOSPHERE 
WITH ASYMMETRIC FORCING 

E. Kalnay and K. C. Mo 


Following a suggestion of Charney (1975, personal communication), Kalnay 
(1977) developed an Iterative scheme to find an unstable stationary solution by 
damping time varying modes. This scheme Is applied to a simple atmospheric 
model In order to find stationary solutions In the presence of asymmetric 
forcing such as orography. 

The model Is a two-level quasi -geostrophic on a P -plane channel with 500 km 
grid size and a domain of 10000 km by 10000 km. It Includes simple parameterl- 
zatlons of long and short wave radiation, surface friction, surface fluxe.s of 
sensible and latent heat similar to those of Charney (1959). The symn«tr^c 
model has no orography and the asymmetric model Is forced by a Gaussian mountain 
located In the center of the channel (Figure 2). 


Figure 1 shows the zonal velocities at the upper and lower levels of the 
stationary solution and the time averaged solution of the symmetric and asymmetrl 
models. Figure 2 shows the upper and lower level streamfunctlons of the 
stationary solution and the t1n« averaged solution of the asymmetric model. 


A linear stability analysis of the four basic states has been carried out 
by successive time Integration of the linearized equations (Brown, 1969) until 
the solution can be expressed as 


')'(t) = 'I'(O) e 


(*^r+1 ®l)t 


Table 1 lists the fastest growing modes and frequencies all four cases. 


Although the time averaged mean fields for both symmetric and asymmetric 
model are far from equilibrium; they are barocllnically much more stable. 
Orographic forcing stabilizes the solutions, especially the time averaged one. 
Wave number 2 Is the fastest growing mode In the time averaged and symmetric 
stationary solutions. In the orographically forced stationary solution, however, 
the most unstable mode Is a combination of wave number 1, 2 and 3, which grow 
at the same rate due to the Interaction between the mountain and the mean flow. 
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Case 

Model 

Solution 

Fastest Growing 
Mode 

Growth 
Rate ■. 
Op (days'*) 

Frequency , 
o^/2tt (day“^) 

1 

symmetric 

stationary 

2 

0.22 

0.036 

2 

symmetric 

time avg. 

2 

0,10 

0,083 

3 

asymmetric 

stationary 

1.2,3 

0.20 

0,062 

4 

asymmetric 

time avg. 

2 

0.041 

0,13 


FIGURE CAPTIONS 


Figure la. Upper level zonal velocities Ui(m/s) for the four cases listed In 
Table 1. 

Figure lb. Lower level zonal velocities U 3 (tn/s) for the four cases listed In 
Table 1. 

Figure 2. The streamfunctlons for the asymmetric model In units of lO^/m^/s. 

(a) upper level, stationary solution, (b) lower level, stationary 
solution, (c) upper level, time averaged solution, (d) lower level, 
time averaged solution. 
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BAROTROPIC INSTABILITY OF WEAKLY NON-PARALLEL ZONAL FLOWS 
L.-O, Merkine* and R. C. Balgovind 


Barotropic instability of weakly non-parallel zonal flows with localized 
intense shear regions is investigated numerically. The numerical integrations 
of the linear stability problem reveal the existence of unstable localized wave 
packets whose spatial structure and eigenfrequencies depend on two parameters 
which measure the degree of supercriticality and the zonal length-scale of the 
shear region. The results Indicate that an unstable global mode solution is 
determined by conditions that ensure the decay of the wave packet at infinity 
and the transition from long to short waves across a turning point (critical 
layer) region which is controlled by non-parallel effects. The controlling 
influence exerted by the weak non-parallel effects on the evolution of the 
instability underlines the weakness of the parallel flow assumption which can 
be used locally, away from critical layers, as a diagnostic tool only. 

The experiments are classified according to the nature of the basic flow, 
namely an intense westerly jet at the center of the channel (Figure 1) or a 
blocking configuration with a split westerly jet (Figure 2). 

The westerly jet results showed that regardless of the type of initial 
conditions (random distribution or a concentrated pulse) instability developed 
whenever the maximum shear was in excess of the critical shear predicted by the 
parallel theory. It assumed the shape of a non-propagating and localized 
envelope modulating propagating waves. The wave-length of these waves was 
location dependent but the frequency of oscillation and growth rate were the 
same for all locations. It follows that the numerical result revealed the 
existence of global modes characterized by a single complex frequency which is 
the eigenvalue of the problem. Figure 3 depicts a structure of the instability 
as a function of x at the center of the channel. The complete spatial structure 
of the unstable wavepack is shown in Figure 4. 

The results of the numerical integrations for a split westerly jet revealed 
again the existence of unstable global mode solutions. The modes describe 
propagating waves modulated by a localized stationary envelope and having a 
single complex frequency. Figure 5 depicts a structure of the instability as 
a function of x at the center of the channel. In contrast to the corresponding 
westerly jet case the deep downstream penetration is absent with the instability 
essentially confined to the supercritical region. The rapid spatial decay of 
the wavepacket is a consequence of the large growth rate. The complete spatial 
structure of the unstable wavepack is shown in Figure 6. 


CONCLUSIONS 

There exists one fundamental difference in the dynamics of instabilities that 
develop in parallel and in non-parallel flows. The parallel flow stability problem 
is invariant to zonal translation since the available energy Is distributed evenly 


* Dept, of Mathematics, Technion - Israel Institute of Technology, Haifa 32000, 
Israel. 
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in the zonal direction. Hence, instabilities can develop without preference at 
any streanwise location. The non-parallel flow stability problem is not 
invariant to zonal translation since the zonal distribution of the available 
energy is not uniform. Consequently, instabilities can develop only at certain 
streamwise locations which are fixed relative to the distribution of the forcing 
functions that drive the basic flow (Merkina, 1982), 

The numerical integration revealed the existence of localized unstable 
global mode solutions whose spatial and temporal characteristics depend on two 
parameters which measure the degree of supercriticality and the length scale of 
the zonal variation of the basic flow. Although we considered a particular basic 
flow, albert with certain appealing features, the basic phenomenon and the 
manner by which it is controlled by these two parameters should be characteristic 
of a wide class of flows possessing localized intense shear regions. 

The results indicate that the complex eigenfrequency of the global mode is 
determined by conditions that ensure the decay of the wavepacket at infinity 
and the transition from long to short waves across a turning point region. The 
dynamics of this transition region which can also be identified as a critical 
layer of the local parallel problem is controlled by the advection of th wave's 
vorticity by the weak meridional velocity of the basic state. The local 
structure of the wave-packet away from critical layers can be determined by 
parallel flow considerations applied locally provided the eigenfrequency is 
known. This underlines the limited usefulness of the parallel theory and 
demonstrates the importance of weak non-parallel effects in controlling the 
dynamics of the instability. 
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Flqure 1. The basic flow for a case when vorticity source is 
neqati ve. 
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Fiqure 2. The basic flow for a case when vorticity source is 
positive. 
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PERSISTENT ANOMALIES OF THE SOUTHERN HEMISPHERE CIRCULATION 

K. C. Mo and J. Shukla 


We have examined the persistent features of large anomalies and determined 
the locations and structures of blocking in the Southern Hemisphere. The data 
set used here contains daily maps of 500 mb geopotential heights for 100 months 
(June 1, 1972 to Nov. 30, 1980), covering from lOS to 90S. The seasonal cycle 
was defined as an 8-year mean and the 8th (annual) and 16th (semiannual) Fourier 
components of the time series at each grid point. Anomalies were defined as 
the difference between the total field and the seasonal cycle for each point. 

The primary conclusions of this study are: 

1. Large anomalies are less persistent in the Southern Hemisphere than in 
the Northern Hemisphere. 

We have calculated the frequency of occurrence of anomalies (^ 150 m or 
< -150 m) which persisted for a certain number of days. There is no evidence 
of any discontinuity in either positive or negative anomaly curves. The peak 
near 2-3 days noted by Charney, Shukla and Mo (1980) due to synoptic scale 
travelling disturbances does not exist here. This indicates that most of the 
large anomalies do not persist more than one day at any grid point. 

2. The local one day lag autocorrelation decay much faster in the Southern 
Hemisphere. Except in the subtropics and Australia, the red noise model fits 
the autocorrelation decay well. 

3. The frequency of occurrence for positive anomalies satisfying the 
criterion (150 m, 6 days) is similar to the map of 8-64 days low pass filtered 
variances, and the frequency distribution for anomalies satisfying the criterion 
(150 m, 1-3 days) represents storm track patterns. 

Figures la and lb show the frequency of occurrence of persistent positive 
and negative anomalies satisfying the criterion (150 m, 6 days) using 100 months 
of data, respectively. Figure la is similar to that presented by Swanson and 
Trenberth (1982) using the criterion (100 m, 5 days). There are three maxima 
located in the southwest Pacific along 55S near 170W and between 140E-150E, 
southeast of South America (65S, 80W) and (503, 80E) Indian Ocean. 

4. The maxima of frequency of occurrence of large positive anomalies and 
large negative anomalies do not occur in the same region. 

Comparison of Figures la and lb shows that the maximum in Pacific (55S, 160W) 
remains in both cases. The one located in the south of southern America shifts 
to (66S, 45W). The maximum at the Indian Ocean does not exist in Figure lb. 

Blocking is defined as an anomaly which satisfies criterion {+ 150 m, 6 days) 
and associates with a pronounced ridge. According to this definition, we have 
identified the blocking events, ridge positions and structures of each event. 

There are 13 events located in the southwest Pacific, four events in the south 
of South America. There is not blocking at (50S, 80E), where the maximum of 
the frequency of occurrence of positive anomalies occurs. However, further 
south at 60S, there are three events in the Indian Ocean. 
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Y. C. Sud and Y. G. Serafini 


This work was aimed to develop a technique to establish the onset of an 
aqricultural drought based on the time evolution of soil moisture predicted 
from the ground water balance. We used a criterion invoking a limiting value 
of 2,5< of the maximum available soil-moisture (Field Capacity - Wilting Value) 
to earmark the onset of an agricultural drought. The soil moisture tendency 
equation requires precipitation and land surface evapotranspi ration as basic 
inputs. The former must be known a priori whereas the latter is introduced by 
suitable parameterizattons for the potential evapotranspi rati on and the moisture 
availability parameter (Mintz and Serafini, 1982), The integrated form of the 
soil moisture tendency equation is found to be; 

w(t) = w^ + 1/X loqe [{ 1 + R exp(-Xw^)} exp(EpXt/CR) - exp(Aw^/R)] (1) 

where X = C = 1 - exp (-«), and R = Ep/[C(P+I) - Ep], 

Here, the time dependent soil moisture w(t) is a function of the potential 
evapotranspi ration Ep, the curvature parameter of the moisture availability 
formulae, the initial soil moisture Wq, the maximum available soil moisture 
w* and the mean moisture source terms; precipitation rate P and irrigation 
rate I. 

By assuming a total absence of all the moisture source terms in an agri- 
cultural drought. Eg. (1) simplifies to yield a time-scale defined as the time 
taken to reduce the mean monthly climatological soil moisture, Wg to the chosen 
wilting value, Wp. Its form is: 

T c/XEp lo% {exp(XWc)-l/exp(XWp)-l} (2) 

The solution of Eq. (1) can be used to predict the time evolution of soil 
moisture to monitor a real drought situation. Furthermore, the constants a, 
w* and Wp can be suitably chosen to best represent the precul iarities of 
the soil type and the vegetation (crops) in a region. 

Alternatively Eq. (1) can be solved to determine the irrigation needs of 
a region. This is equal to the difference of the calculated monthly moisture 
source term (required to maintain a certain mean available soil moisture) and 
the precipitation. 

Using the constructed fields of global mean '’'onthly climatological values 
of soil moisture (Mintz and Serafini, 1982), the time period for the inception 
of an aqricultural drought was obtained from Er^. (2) for each month. Sample 
plots for the months of January and July are shown in Figs. 1 and 2. The 
corresponding climatological soil moisture values are shown in Figs. 3 and 4, 
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DEPENDENCE OF PRECIPITATION ON EVAPOTRANSPIRATION 
IN THE CONTINENTAL UNITED STATES 

X, C. Sud and M* Fennessy 


The Influence of soil mosture and Its dependent evapotranspiratlon on mean 
inonthly precipitation was examined by analyzing forty years of observed precipi- 
tation and surface temperature data. The data comprised of actual observations 
averaged over each month at 261 local climatic data (LCD) stations distributed 
throughout the continental United States covering the period from 1939-1980, 

The actual evapotranspiratlon was obtained by first calculating the mean 
monthly potential evapotranspiratlon using empirical functions of Thornthwalte 
(1948) and Wlllmott (1977) and then multiplying It by the soil moisture dependent 
moisture availability parameter (Nappo, 1975). This evapotranspiratlon In con- 
junction with the observed precipitation provided the net surface water balance 
In the soil moisture tendency equation. The soil moisture tendency equation 
was Integrated to yield the soil moisture and evapotranspiratlon at the 
beginning of each month. 

The observed mean monthly precipitation was correlated with the antecedent 
fields of a) observed mean monthly precipitation, b) soil moisture at the be- 
ginning of the month, and c) evapotranspiratlon at the beginning of the month. 

Statistically significant correlations were found In the drought prone 
western Freat Plains region In September and October. 

Fig. 1 shows the correlation coefficients between precipitation and antece- 
dent soil moisture for the month of October. The Great Plains regions was then 
further analyzed by constructing three sets of contingency tables between the 
correlated fields. These statistics reaffirmed the existence of significant 
positive correlations between precipitation and other fields. A contingency 
table between soil moisture and precipitation Is shown In Table 1. The strength 
of diagonal from top left to bottom right represents the positive correlation. 
CHI-square (v^) values at the bottom of the table reflect the strength of the 
signal (s^ » for purely random distribution). 

Based on this study we conclude that In the drought prone western Great 
Plains region a statistically significant relationship exists between antecedent 
soil moisture and the precipitation, particularly In the later part of summer. 
This suggests the Importance of soil moisture In maintaining rainfall anomalies. 
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Normalized Contingency Table Between Antecedent Soil Moisture and Precipitation 
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OCTOBER NCC LCD COR. CO. PRECIP- ~ SOIL MOISTURE 



Fig. 1. Large print is for correlation coefficients having more than 
801 significance for a 2-tailed T-test. 
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DESIGN OF AN INTERACTIVE BIOS.'HERE 
FOR THE GLAS GENERAL CIRCULATION MODEL 

Y. Mintz, P. J. Sellers and C. J. Wlllmott 


Research has begun on an Interactive biosphere for the GLAS general 
circulation model, so that the calculations of the land-surface evapotrans- 
piration and sensible heat flux will be more realistic. This is important 
because sensitivity experiments with general circulation models have shown 
that the simulated fields of rainfall, temperature and motion of the 
atmosphere are highly sensitive to the transfers of latent and sensible 
heat at the land surface. 

Water and energy transfers at the land surface depend on vegetation 
morphology and physiology. In the model biosphere, the vegetation morphology 
will control these transfers through the leaf area density as a function 
of height, Ln(Z), and through the root length density as a function of 
depth, RTo(Z), as expressed in a discretized form for each of the model 
grid areas. Lo(Z) will influence (a) the aerodynamic resistance to the 
latent and sensible heat transfers; (b) the transmission, jbsorption and 
emission of radiation energy by the canopy and underlying ground; and (c) 
the interception, evaporation and throughfall of rain and snow. The vege- 
tation physiology will control the transfers of water and energy through 
the stomatal resistance, the xylem resistance of the stems, and the root 
cortex resistance. The stomatal nssistance will be a function of the leaf 
water potential, short-wave radiation intensity, leaf temperature, and 
humidity of the air. 

In the initial formulation, Lo(Z) and RTo(Z) are being prescribed as 
functions of latitude and longitude and the season of the year, as known 
from ecological observations in the various vegetation formations. Later, 
the phenological changes of Ld(Z) and RTq(Z) in the deciduous forests 
and grasslands will be made interactive with the atmospheric variables 
and the soil moisture as calculated by the GCM. Finally, the vegetation 
formations themselves will be made interactive with the atmospheric and 
soil moisture conditions. 

The realism of the model biosphere will be evaluated (1) through 
short period one-dimensional comparisons with measured values of the local 
atmospheric forcing and measured values of the vertical fluxes of water 
and energy; (2) through long period comparisons of simulated and observed 
catchment water budgets; and (3) through comparisons of the simulated and 
observed surface temperatures, surface albedos, snow cover, live biomass, 
and the water balances of the large river basins, when the model biosphere 
is forced with atmospheric variables taken from the FGGE level 3-B and 
2-C data sets. 

The accompanying figure illustrates (for the case of a rainforest) the 
procedure whereby the morphological and physiological characteristics of the 
vegetation are transformed into the parameters which govern the transfer of 
water from the root zone of the soil to the atmosphere. (For further 
details, see NASA Technical Memorandum XXXXX, January 1983.) 





153 


PRECEDING PAGE BLANK NOT FILMED 







A. Tropical rainforest vegetation; continuous canopy of tall 
trees above continuous under-storey* 

B. Leaf area density S , stem area density ES3 , and 
root length density CQIiT] , as functions of height for the 
tropical rainforest shovm in Figure A. 

C. Electrical resistance analogue of the tropical ralnforesr 

represented in Figure B» F raintorest 



INFLUENCE OF EVAPORATION ANOMALIES IN SEMI-ARID REGIONS ON JULY CIRCULATION 

Y. C. Sud and M. Fennessy 


In this study, two 47-day siinulatlons were made with the general circulation 
model (GCM) of the Goddard Laboratory for Atmospheric Sciences (GLAS) starting 
from the observed initial conditions for June 15 (OOZ) 1979 as produced by NMC 
analysis. In the first simulation, the soil moisture was normally initialized 
and then allowed to vary by the model generated surface water balance at each 
grid point. In the second simulation, the soil moisture was maintained at 
zero, yielding NO evaporation in each of the four chosen semi -arid regions 
shown in Table TT Elsewhere, the soil moisture was calculated as in the first 
simulation. The two runs were identical in all other respects. 

The mean July precipitation in three out of four regions increased. These 
results are summarized in Table II, In the monsoonal regions of the Northern 
Hemisphere, particularly the Indian subcontinent, these results accord with a 
recent study by Shukla and Mintz (1982). 

In each region the analysis revealed an increase in the mean monthly 
diabatic heating, vertical velocity and moisture convergence fields. Some 
increase in the moisture convergence due to an increase in the sensible heat 
flux (in the absence of evaporation) was to be expected. However, the overall 
influence of an overcompensation of the moisture deficit is truly surprising. 

This study has supplemented an earlier study by Sud and Fennessy (1982) 
wherein the albedo was increased in the very same regions. The two studies 
together lead to the following conclusion. Of the two surface forcing parameters, 
the albedo and the soil moisture, when occurring in small semi-arid regions 
bordering deserts it is the influence of the albedo which tends to reduce the 
rainfall by producing diabatic cooling and subsidence. The evaporation anomaly 
influence, on the contrary, may largely compensate the moisture deficit by 
inducing moisture convergence caused by a corresponding increase in the sensible 
heat flux. It is, therefore, apparent that Charney's hypothesis which stipulated 
albedo increase as the primary mechanism for the feedback of a desert on itself 
is provided further support by this study. The role of soil moisture, though 
important, does not have such a positive feedback effect. 

Charney et al . (1977), in their study with prescribed soil moisture, found 
a much larger influence of albedo changes on the precipitation as compared to 
Sud and Fennessy (1982). It can be argued that a part of the reason for the 
discrepancy between the two results is attributable to the prognostic soil 
moisture calculation in the latter study, which presumably had a negative feed- 
back effect on the ensuing precipitation. 

In interpreting the guantitative significance of these results, the limita- 
tions of the GCM must be recognized, particularly in parameterizations of the 
planetary boundary layer and convective clouds. 


155 



REFERENCES 


Charney» J, G., W. J. Quirk* S. H. Chow, and J. Kornfield, 1977: A comparative 
study of the effects of albedo change on drought In semi -arid regions. ^ 
Atmos. Scl. , 34 , 1366. 

Shukla, J., and Y. Mintz, 1982: The Influence of land surface evapotranspiratlon 
on earth's climate. Science , 215 , 1498-1501, 

Sud, Y. C., and M. Fennessy, 1982: A study of the Influence of surface albedo 
on July circulation In semi -arid regions using the 6LAS 6CM. J. of Clima- 
tology . 2, 105-125. 


Table I 


original page is 

OF POOR QUALITY 


Areal Pctent of Wo Evaporation Anoaialy Wegiona 


Region 

Nueber of 
Grid Points 

Latitude 

Longitude 

Sahel 

22 

12*M-16«N 

17.5*H-47.5*B 



16*H-20*N 

17.5«W-37.5*B 

Thar Desert 

4 

24*M-32»N 

67.5*E-77.5*B 

Border 




Brasil 

10 

4*S- 8*S 

32.5*W-^37.5*W 



4*S-20*S 

37.5*H-42.5*H 



4*S-24*S 

42.5*H-47.5*H 

Has tern Great 

S 

32*N-4B*N 

97.5*H-107.5*M 

Plains 










Rsqlonal Diff«rqne«i {AnoMl 



ORtGINAL PAGE IS 
OF POOR gUALflY 



Mot«: percentage changes are based on the siaulated values in the control 

























ORlGSi^JAL PASS IS 
OF POOR QUALITY. 


AN OBSERVATIONAL-DATA BASED EVAPOTRANSPIRATION FUNCTION 
FOR A GENERAL CIRCULATION MODEL 

Y. C. Sud and M. J. Fennessy 


As described by Carson (1981), the evaporation calculation in various 6CMs 
is of the form 

E(Tg) B (w,w*) Ep(Tg) , (1) 


where E is the actual evaporation, and Ep is the potential evaporation from the 
surface of the ground at temperature Tg and B(w,w*) [hereafter 3] is the evapo- 
ration coefficient and w,w* are actual and saturation field capacities. 

On the other hand, field measurements [e.g. Davies and Allen, 1973] usually 
give the moisture availability coefficient M(w,w*) [hereafter M] which is 
defined as: 

M = E(Tg)/Ep(Tg*). (2) 

Here Tg* is the temperature of the saturated ground interacting with the atmo- 
sphere above. In order to use (1) to obtain E(Tg), a 3 function based on 
observational data, such as available in the form of tfje M function, is needed. 
Thus far, simple linear 3 functions have been used in various GCMs without even 
properly distinguishing between M and 3. 

For the experiment (Davies and Allen, 1973) a single relationship, M = M 
(w,w*), was found by Nappo (1976) by averaging the measurements. A similar 
approach has been followed to derive the 3 function. A well known version of 
the GLAS GCM was modified to obtain two ground temperatures and two potential 
evaporations, from the naturally dry and saturated ground respectively. The 3 
function values were derived from thj daily mean evaporation fields (1). One 
hundred ensembles of egually spaced soil moisture between a minimum, of zero 
(dry), and a maximum of one (saturated), were used to calculate 3 from the last 
31 days of simulation. 

An examination of the 3 curve led us to a functional form of the type: 

3 = bWf + cWf^ + aw^^/ (Wjr^ + e). (3) 

The constant e was introduced to produce a small curvature to the desired curve 
near the low soil moisture range of Wf 0. Values of the constants are given 
in Table 1. The corresponding 3 curve is shown in Fig. 1. 

It must be noted that the 3 function derived in this analysis still gives 
a simple evapotranspi ration function over global lands. 
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a 

b 

c 

Figure 2 
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.236 

o 

t-H 

• 

1 

Figure 3 

0.100 

1.09 

.087 
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lA STRATOSPHERIC VERSION OF THE 6LAS CLIMATE MODEL 
J. Abeles and D. Randall 


The stratosphere is a region of great Importance to the earth's climate 
system. In the stratosphere, the photochemical ly and radiatively active sub- 
stance of ozone provides a shield which is vital to life on the planet. In 
addition, the large-scale dynamical processes occurring in the stratosphere can 
have a strong impact on the tropospheric circulation. The mere presence of a 
stratosphere is helpful in correctly modeling the troposphere, because the 
removal of the model top to a very high altitude where significant dissipation 
processes come into play mitigates false wave reflections due to the rigid top 
boundary condition (Lindzen et al., 1968; Kirkwood and Oerome, 1977; Mechoso et 
al. . 1982). 

To successfully model the troposphere-stratosphere system, it is necessary 
to take full account of the tropospheric processes that significantly influence 
the stratosphere. These involve a very wide range of space and time scales. 

They include vertically propagating planetary waves, the influence of tropospheric 
clouds on the upwelling long-wave radiation at the tropopause, deep moisture 
transport by the Hadley cell, and the possible forcing of the semidiurnal tide 
in the lower stratosphere by cumulus heating. All of these processes are 
strongly Influenced by the release of latent heat, which is determined by the 
parameterized processes of the planetary boundary layer (PBL) and cumulus convec- 
tion. 


We are conducting a study to examine the impact of these processes on the 
stratosphere by directly comparing the tropospheric forcing of the stratosphere 
and the simulated stratospheric circulations produced by two models whose tropo- 
spheres differ only in their treatment of diabatic heating. Although existing 
observational data on the stratosphere are too meager to allow a detailed 
comparison of model results with observations, particularly for the tropics and 
for seasonal variations (Hudson and Reed, 1979), gross comparisons are feasible 
and are an important objective of our study. 

The first model is based on the potential -enstrophy-conserving GLAS Climate 
Model (hereafter the GLAS model). This model incorporates the tropospheric 
physical parameterizations of the model described by Shukla et al . (1982), 
including a cumulus parameterization based on the theory of Arakawa (1969), and 
a planetary boundary layer (PBL) parameterization based on the work of Katayama 
(described by Arakawa, 1972). In the past, the GLAS Climate Model has used nine 
layers, of equal thickness in sigma, with the model's top at 10 mb. Although 
the 10 mb surface is well within the stratosphere, the top model layer was so 
deep (on the order of 15 km) that the model could not resolve any stratospheric 
structure; we effectively had a tropospheric model. 

During the past year, we have created a new, stratospheric version of the 
model. This new model differs from the old in the following respects: 

1) the model's top is at 1 mb; 

2) there are 15 layers, in all; 
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3) a pressure coordinate system is used above 100 mb; 

4) the eight sigma-layers below 100 mb are of equal thickness, to accommo- 
date certain restrictive assumptions built into the GLAS cumulus para- 
meterization, but the seven pressure-layers above 100 mb can have 
arbitrary thicknesses; 

5) ozone has been introduced as a prognostic variable, and the simple 
ozone photochemistry parameterization of Schlesinger and Mintz (1979) 
has been incorporated into the model; 

6) the Wu-Kaplan long-wave radiation parameterization has been replaced by 
the Schlesinger-Mintz parameterization; 

7) the solar radiation parameterization has been generalized by R. Davies, 
of Purdue University, to allow a general distribution of ozone, and an 
arbitrary number of model layers with arbitrary thicknesses. 


The second model is based on the UCLA model developed by Arakawa and col- 
laborators at UCLA, and currently running at GSFC (Abeles, Corsetti, and Randall, 
this volume). This model incorporates the Arakawa-Schubert-Lord-Chao cumulus 
parameterization and the Randall PBL parameterization. 

The two models are essentially identical except for their differing para- 
meterizations of cumulus convection and the PBL. This will allow us to study 
the effects of these two key tropospheric physical parameterizations on strato- 
spheric circulations, by performing numerical experiments in which the two 
models are run with identical initial and boundary conditions. We will place 
special emphasis on comparison of the model results with respect to vertically 
propagating Planetary waves, upwelling long-wave radiation at the tropopause, 
deep moisture transport by the Hadley cell, and cumulus forcing of the semidiurnal 
tide of the lower stratosphere. 

Currently, the stratospheric GLAS model has been debugged, and a January 
simulation is under way. A January simulation has been completed with the 
15-level UCLA GCM. Comparison of the results obtained with the two models will 
begin in the near future. 
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SEASONAL CYCLE SIMULATION 


L. Marx and D. Randall 


The GLAS climate model Is being used In studies of the predictability of 
monthly and seasonal means » and In tests of the sensitivity of the general 
circulation to boundary forcings such as anomalies of sea surface temperature 
or soil moisture. Recently* the model has been run for over two simulated 
years, to test Its response to nonInteractIve annually varying boundary condi- 
tions. 

The past two years have seen major Improvements In model performance, as 
a result of several revisions of the model's numerical schemes and physical 
parameterizatlons. The principle changes to the model are: 

1. New boundary layer parameterization; 

2. New parameterization of evapotranspi ration from land surfaces; 

3. New parameterization of the Interaction between clouds and radiation; 

4. New scheme to maintain computational stability near the poles; 

5. New vertical differencing scheme; 

6. New tape format; 

7. New postprocessing package. 

These changes are described by Randall (1982). The model code has been completely 
rewritten for Improved speed and readability, and we have produced a comprehensive 
documentation. The model Is running on the new Cyber 205 computer at 6SFC. It 
uses approximately 3 cpu minutes per simulated day. 

A completely new postprocessing package has been designed to aid In the 
study of the model results. The synoptic records described above are written 
to tape In sigma coordinates. A program has been created to read one or more 
of these sigma-history records. Interpolate to constant pressure surfaces, and 
write the results to tape. Another program Is used to create tapes containing 
time-averages of these pressure history tapes. An elaborate plot package has 
been developed to generate plots of Individual pressure history records, time 
averages, or differences between time averages. Finally, grid-point history 
records can be sorted by grid point, and plotted as simple plots against time 
(e.g., the record of precipitation at a grid point) or as time-height contour 
plots (e.g., of total diabatic heating). 

Two changes were made to the boundary conditions for this run. First, the 
topography was revised to assign more realistic heights to Greenland. Second, 
the seasonally varying ground wetness data of Mlntz and Serafinl (1983) was 
Included as an additional boundary condition, I.e., the ground wetness was 
made non-lnteractlve. 
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The first 300 days of the run were carried out on the Cyber 203 at the 
Langley Research Center, and the remainder was completed on the Cyber 205 at 
6SFC. The model Incidentally served as a benchmark for the newly Installed 
Cyber 205. The running speed Improved by a factor of 2.5 on the 
Cyber 205, as compared with the Cyber 203. 

The model was Initialized with the observed state of the atmosphere for 
OOZ, November 15, 1978. The Initial ground temperatures for both the grid box 
and the small saturated plot were assumed to be equal to the surface air temper- 
ature. The Initial snow boundaries were obtained from the climatological data 
of Matson (1977), and the Initial snow depths were specified as functions of 
the surface albedo. 

The results of the run, represented by monthly mean plots of several hundred 
prognostic and diagnostic fields. Indicate realistic responses to the variations 
In the boundary conditions. There are realistic annual variations In the sea 
level pressure, geopotential height, and other fields over the two years of 
Integration. The global and hemispheric means of potential temperature, kinetic 
energy, precipitation and evaporation also show realistic seasonal variations. 
Some results are shown In Figs. 1 and 2. 

Perhaps the most serious deficiencies In the results are the familiar 
"cold-pole" problem, and an excess of low-level stratiform cloudiness. These 
problems are being addressed In ongoing research. 
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AN OBSERVATIONAL STUDY OF LARGE SCALE ATMOSPHERIC ROSSBY WAVES 
R. S. Llndzen, D. M. Straus, and B. Katz 


Analyzed global 500 mb data from the National Meteorological Center for 
the period 1 January 1975 to 28 February 1975 are used to determine phase and 
amplitude as functions of time for westward traveling Rossby waves. The 
methodology consisted of projecting the data (at each time) onto Hough functions, 
followed by band pass filtering to retain only westward propagating fluctuations 
In a broad period range. The first three meridional and zonal wavenumbers were 
studied, and In most cases Rossby waves with phase speeds close to theoretical 
expectations were found. 

It is clear that steady phase propagation does not occur (see Figure 1 as 
an example). Typical Intervals between breakdowns in uniform phase progression 
are comparable to the wave period. Such breakdowns are generally associated 
with minima In wave amplitudes. This picture Is consistent with episodic wave 
generation and subsequent frictional decay. While inviscid, steadily propagating 
Rossby waves are not expected to Interact with the mean flow (averaged over 
Rossby wave periods), the same Is not true for damped waves. Our analysis 
clearly demonstrates the damped character of atmospheric Rossby waves. How 
Important such waves are to atmospheric transport remains to be determined. 
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Figure 1. The behavior of the (1,2) Rossby mode (zon’l wavenumber 1, 
one zero crossing between the poles). The >hase convention 
is that phase Increases for westward propagv^tion. The solid 
line in the lowest box is the amplitude (use the left-hand 
scale for the height field, the right-hand sc vie for the 
wind field). The dotted line in the lowest bc;< is the 
amplitude of the average Hough mode, where the average was 
taken over rolling blocks of length 6.5 days. The middle 
box gives the percentage of the total variance if the Hough 
mode projection over a particular block that was explained 
by the Rossby wave. 



ON THE ROLE OF THE SEASONAL CYCLE 
D. H. Straus 


The periodic change In the general circulation and local climate (i.e. 
the seasonal cycle) stands out as the most noticeable climate signal on time 
scales less than a decade or so. The resulting non-stationarity in the mean is 
not taken into account in the conventional statistical framework in which 
transients are considered to be departures from the time mean. A more natural 
definition of transients is the departure from the seasonal cycle. 

The seasonal cycle can be defined in terms of a few orthogonal polynomials, 
as described in detail by Straus (1983). This is particularly useful when only 
short time series are available, so that ensemble averaging is not possible. 

The choice of polynomials depends on the length of the time series in question; 
the seasonal cycle is the projection of the original series onto the subspace 
of the chosen polynomials. This insures that the time average of any product 
of variables over the record length equals the sum of the time average of the 
product of the seasonal cycle contributions and the time average of the product 
of transient terms, with no cross terms. If the time average is taken over a 
period less than the record length, cross terms (or interaction terms) are 
involved and must be interpreted. 

An interesting example of the importance of interaction terms is presented 
in Straus (1983). The 200 mb eddy momentum flux (zonally averaged) is examined 
from a 7-year record. The harmonic sines and cosines form a natural basis set 
with which to describe the seasonal cycle, which then consists of the 7-year 
time mean and oscillations of periods one year and 1/2 year. The remaining 
fluctuations are then divided up (using Fourier analysis once more) into medium 
frequency fluctuations (periods of 7.2 days down to the Nyquist frequency of two 
days) and low frequency fluctuations (involving periods of 7.2 days to 7 years, 
excluding the annual and semiannual cycle). The focus here is on the interannual 
variability of the eddy momentum flux. Winter (December through February) 
averages of the flux were examined for each of the 7 years. Since the time 
averaging period (3 months) is not identical to the record length over which 
the seasonal cycle is defined (7 years), interaction terms are present. Thus, 
for any one winter, the total eddy momentum flux is the sum of a seasonal 
cycle contribution, a low frequency contribution, a medium frequency contribu- 
tion, and interactions between these different time scales. It turns out that 
the only interaction term of importance is that involving the low frequency 
signal and the seasonal cycle. This term explains most of the interannual 
variability seen in the total eddy momentum flux (see Fig. 1. This indicates 
that the low frequency anomalies themselves are not primarily responsible for 
the interannual variability of the fluxes. Rather, it is their positioning in 
time with respect to the seasonal cycle that is important. It is hoped that 
these concepts will prove useful in the study of mid-latitude blocking and the 
phenomenon of El Nino. 
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Figure 1. Winter averaged, zonal 1y averaged eddy momentum flux 
at 200 mb as a function of year (labelled 1 to 7) at 
different latitudes. The heavy solid line is the 
total flux, the light solid line is the flux due to 
the low frequency components, tie dotted line is the 
flux due to the medium frequency components and the 
dashed line gives the flux due to the low frequency- 
seasonal cycle interactions. When a separate momentum 
flux scale is given as the right side of a panel, it 
refers to the total flux only. 







INFLUENCE OF SOIL MOISTURE ON CLIMATE 
J. Stiukla and Y» Mintz 


The Goddard Laboratory for Atmospheric Sciences (O'.fS) Climate Model has 
been used to conduct a numerical experiment to determi ^i? t'e influence of soil 
moisture on climate (Shukla and Mintz, 1982). The model Wiij run for 60 days 
with two assumed extreme global conditions of soil moisture, dry and saturated; 
for dry soil evaporation from the land surfaces in zero, and for saturated (wet) 
soil the evaporation equals the potential evapotranspi ration. The differences 
in the simulated atmospheric circulations were enormous. For example, over most 
of North America and Europe the simulated rainfall for dry soil was about 40 
percent less compared to the wet soil. The study showed that evaporation from 
land is an important factor in determining the rainfall over land. The results 
from this experiment suggest the important role of soil moisture in the inter- 
annual variability of climate and indicate that global observations of soil 
moisture can provide potential predictability for monthly and seasonal anomalies 
in the atmospheric circulation. 


DESCRIPTION OF FIGURES 

The top panel shows the annual and global mean water budget. The annual 
mean evaporation over land (41 cm) is about 60% of the annual mean precipitation 
over land (72 cm). The runoff from land to ocean is 31 cm. 

The lower three panels show the results of the numerical experiment. 


(a) Ground temperature: Difference [dry soil (no evaporation) - saturated 

soil (maximum evaporation)] for July in ("C). 

° Positive values mean a warmer ground temperature. When land is dry, 
all of the solar radiation goes to heat the ground (and none to 
evaporate the v^ater) and ground temperatures are very high. 

(b) Sea level pressure: Difference [dry soil (no evaporation) - saturated 

soil (maximum evaporation)] for July in (mb). 

° Maximum re^luctions of pressure occur over land for the dry soil case 
in the Northern Hemisphere. A compensating increase in pressure occurs 
over oceans. 

(c) Precipitation: Difference [dry soil (no evaporation) - saturated soil 

(maximum evaporation)] for July in mm/day. 

* In the case of dry soil, there is a significant decrease of precipi- 
tation over land over most of the globe except for the monsoon regions. 
Over the Indian monsoon region, large surface heating creates very 
large moisture flux convergence (moisture comes from the adjoining ocean 
and leads to larger precipitation. 
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Globally averaged precipitation (over land and ocean) in the dry soil 
case is reduced by about 18% of the saturated soil case. Globally 
averaged precipitation over land only is reduced by about 54% of the 
saturated soil case. 


REFERENCE 

Shukla, J.» and Y. Mintz, 1982: Influence of land-surface evapotranspi ration 

on the earth's climate. Science, 215, 1498-1501. 


176 


original page is 

0^ POOR QUALITY 










INFLUENCE OF EQUATORIAL SST ANOMALIES (EL NINO) ON THE CLIMATE OF NORTH AMERICA 


J. Shukla and M. Wallace 


A number of recent observational and theoretical studies indicate that the 
equatorial Pacific SST anomalies, also referred to as El Nino events, can 
produce significant climate anomalies over North America. This phenomenon is 
illustrated schematically in the bottom panel where the shaded area represents 
a zone of enhanced precipitation associated with warm SST anomalies which force 
a wave pattern along a great circle path. 

A numerical experiment was conducted (Shukla and Wallace, 1982) with the 
GLAS climate model to test the validity of the above hypothesis with January 
initial conditions based on observed data and an equatorial Pacific SST anomaly 
based on the recent analysis of Rasmusson and Carpenter (1982) shown in the top 
panel. The middle panel shows the difference between the anomaly and control 
integrations averaged for days 11-30 for geopotential height at 300 mb. Positive 
and negative values of height differences Indicate areas of warmer and colder 
air temperatures, respectively. It is seen that the maximum changes occur over 
North America. These results further confirm the earlier suggestions based on 
observational and simple model studies, and highlight the Importance of equatorial 
Pacific SST anomalies in variability and predictability of monthly and seasonal 
mean atmospheric circulation over North America. 


DESCRIPTION OF FIGURES 

Top Panel: Observed SST anomalies averaged for Nov.* Dec., and Jan., for four 

El Nino episodes (Rasmusson and Carpenter). 

Middle Panel: Difference (anomaly-control) of 300 mb height averaged for days 

11-30. The thick solid contour denotes the zero value. 

Bottom Panel: Schematic representation of upper tropospheric teleconnection 

patterns associated with tropical rainfall anomaly over the shaded region. 
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QUASI-STATIONARY STATES AND EDDY-MEAN FLOW INTERACTION IN A WIND-FORCED 

BAROCLINIC OCEAN MODEL 

K. C. Mo, E. Kalnay and D. E. Harrison 


We present results from a square basin, steadily forced non-linear two- 
level quasi geostrophic ocean model experiment, designed to study the effects of 
mesoscale motions on the mean flow. A conventional calculation of the behavior 
resulting from spin up from rest leads to the appearance of baroclinic eddies 
after about 800 days of integration, followed by a statistical equilibrium in 
which the eddies are prominent. The mean statistical equilibrium flow (SEQ) is 
similar to others which have been discussed previously, although we employ a 
different numerical method. We also present two quasistationary (QS) solutions 
to the same model ocean system, found using an iterative scheme that damps out 
time dependent modes, thus suppressing the appearance of eddies. The QS 
solutions result from convergence of the eddy-suppressing scheme started from 

different initial conditions. We find that the RMS value over the basin of the 

time derivative of the potential vorticity (R) of the QS solutions is: 1) 

locally minimum in the sense that small perturbations to the QS solution increase 
R and 2) Rqs is much smaller than the RMS of the other terms in the model 
potential vorticity equation. We find also that Rqs « RsEQ* 

When the first guess for the iterative scheme is a state of rest, the 

system converges to a quasistationary solution denoted QSl. From any other of 
our initial guess states, the scheme converges to a second quasi stationary 
solution, QS2. QSl resembles the state of wind driven ocean attained after 
about 350 days integration from rest, whereas QS2 resembles the state of the 
ocean after about 650 days of integration. 

When QSl is used as the initial condition in a conventional time integration, 
a recirculating gyre mode grows linearly with time for about 400 days, even 
after the mode has attained an amplitude of the same order as QSl. This linear 
growth ceases only when the flow becomes baroclinically unstable and eddies 
start to develop. Thereafter the SEQ solution is recovered. We discuss the 
dynamics of this resonant, linearly growing mode; it is not simply an inertial 
gyre. QS2 resembles QSl plus the linearly growing mode after 400 days. When 
QS2 is used as initial condition, linear growth also occurs but the flow becomes 
baroclinical ly unstable much faster. 

The role of the baroclinic eddies in this model ocean is studied in terms 
of the differences between SEQ, the statistical equilibrium solution obtained 
in the presence of eddies, and QS2, the second and most intense quasistationary 
solution. Overall, the eddies weaken the mean upper layer flow in the recircu- 
lation, and alter the mean lower layer circulation signficantly. 


FIGURE CAPTIONS 


Figure A: 

Figure B: 

Figure C; 
Figure D: 
Figure E: 
Figure F: 


The upper level streamfunctlons 4»i for SEQ In 100 m^/s; contour 
Interval 4000 mvs. 

The lower level streamfunctlons i |»3 for SEQ In IOC m^/s; contour 
Interval 800 m^/s. 

Same as Figure 6a but for QSl. 

Same as Figure 6b but for QSl, 

Same as Figure 6a but for QS2. 

Same as Figure 6b but for QS2. 
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IV. CYBER CONVERSION 


EXPERIMENTS WITH THE UCLA GCM 
J. Abeles, T. Corsettl, and D. Randall 


1. MOTIVATION 

The UCLA General Circulation Model, developed over a period of two decades 
by A. Arakawa and collaborators at UCLA, has many features in common with the 
current GLAS Climate Model, yet it offers a number of significant advantages. 

Most importantly, it includes more realistic physical parameterizations, parti- 
cularly for the planetary boundary layer and cumulus convection. It also employs 
improved finite-difference schemes for horizontal advection and geostrophic 
adjustment. It includes the mixing ratio of ozone as a prognostic variable, 
whose distribution is modified by advection and a simple photochemistry parameter-^ 
ization. Finally, the program structure is flexible as to horizontal and vertical 
resolution. The model description is summarized in Table 1. We are in the 
process of performing simulations with the model, to evaluate its ability to 
simulate the general circulation. 


2. NUTS AND BOLTS 


An effort to optimize the model on the Cyber 205 is ongoing. The original 
code was taken from Monterey, where limited vectorization had been accomplished 
on a Cyber 203, The work began on Cybernet (a commercial network) and has 
c„,tinued on the NACF 205. Numerous problems were encountered during the 
conversion effort, including compiler bugs detected during the execution of the 
model . 

The initial goals included creating "logical unit" subroutines, making the 
code more readable, and creating subroutines to be shared between the physics 
and hydrodynarfiics , as well as top-down vectorization of the hydrodynamics. 

Table 2 shows the spied tncre^ise to date by vectorization combined with scalar 
optimization of various reutirigs. Vectorization of the physics has also begun. 

The model history data is interpolated to pressure surfaces, converted 
from full precision Cyber words (64 bits) to full precision Amdahl words (32 
bits), and trarisferred to the front-eitd for further analysis. A program to 
process grid point history data has been completed; as a result, we will be 
able to put out a very detailed record of each simulation, for selected grid 
points. 


3. RESULTS 


The initial problems have been oyercome, and two model simulations have 
begun. A i5«level model with its top at i mb has completed a January simula- 
liffi, and a 9-1evel model with its top at BO mb has completed a simulation which 
fXt^^ndS from December 31 to the end of July. Figs. 1-4 show the simulated 
200 mb geopotential height, PBL depth, total cloudiness, and PBL stratus inci- 
d@nee, for July of the 9-level simulation. The 2j00 mb height field shows 
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strong, very realistic planetary wave activity. The PBL depth shows maxima 
over the warm coastal currents of the Northern Hemisphere, and minima of land 
where the nocturnal values are small. The total cloudiness has a realistic 
geographical distribution, and the globally averaged value is also realistic. 
The PBL stratus incidence shows subtropical maxima, in agreement with observe 
tions. 
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Table 1 

MODEL DESCRIPTION 


A. Prognostic variables: surface pressure, planetary boundary layer depth, 
zonal wind, meridional wind, potential temperature, specific humidity, 
ozone mixing ratio, ground temperature, snow depth. 

B. Vertical structure; 

1. The stratosphere uses a pressure coordinate system 

2. The troposphere uses a sigma coordinate system 

3. The lowest model layer is the planetary boundary layer 

4. Arakawa-Suarez-Lamb vertical differencing. 

C. The horizontal and vertical resolution are fully variable. 

D. Horizontal differencing: Fourth-order potential -enstrophy-conserving 

momentum advection scheme developed by Kenji Takano, and fourth-order 
advection of potential temperature, moisture, and ozone. 

E. Leapfrog-Matsuno time differencing. 

F. Physical parameterizations: 

1. Arakawa-Schubert cumulus parameterization 

2. Randall planetary boundary layer parameterization 

3. Schlesinger-Katayama radiation parameterization. 
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Table 2 

15-LEVEL MODEL TIMINGS (minutes per simulated day) 


INITIAL NOW 


Machine 

203 

205 

205 

205 

Compiler 

1.5 

1.5 

1.5 

2.0 

Compiler Options 

o.u.v 

o.v 

NONE 

MIXED 

Hydrodynamics 

15.6 

17.2 

55.2 

4.0 

Physics 

11.5 

9.5 

13.2 

7.0 

TOTAL 

27.8 

27.3 


11.4 
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VECTOR I ZAT I ON OF THE PHYSICS IN THE GWS FOURTH ORDER MODEL 

Lawrence L. Takacs 


In order to utilize the high-speed vector processing capability of the 
Cyber 205 computer, the 'physics' of the GWS fourth order general circulation 
model has been greatly restructured and modularized. The physics of the model 
consists of the determination of the ground condition, the calculation of in- 
coming short-wave and outgoing long-wave radiation, the effect of cumulus 
convection, the determination of precipation, etc,, and was previously calculated 
one grid-point at a time with vertical integration at each grid-point where 
appropriate. The new version of the code streams through all longitudinal 
grid-points at the same time and, where appropriate, the vertical index as well. 
This restructuring and vectorization of the physics has produced an overall 
increase in speed by a factor of 7.6 over the Cyber scalar version when long- 
wave radiation is called, and by a factor of 3.8 otherwise. 

The main body of the physics code is located in subroutine C0MP3. The 
cumulus parameterization has been moved to a new subroutine called CUMULO, 
while the calls to the shortwave and longwave radiation routines are done in 
the new routine RADIO. The vectorization of the code is an ongoing project, 
and further improvements and enhancements of the model are currently being 
investigated. 
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VECTORIZATION OF THE HYDRODYNAMICS IN THE 6LAS FOURTH-ORDER MODEL 

Ramesh C. Balgovind 


The Hydrodynamics of the GLAS Fourth-Order General Circulation Model 
FORTRAN code has been converted from scalar to a vector form. This in order 
to utilize the high-speed vector processing capabilities of the Cyber 205 
computer. The Hydrodynamics part of the model consists of the advection» 
Fourier filterings, Shapiro filtering and the calculation of pressure to the 
kappa power. 

In the conversion effort our major objectives were to; 

1) vectorize as many DO loops as possible, 

2) minimize division and the use of Q8 calls, 

3) avoid descriptors whenever possible, 

4) keep the variable resolution, 

5) form link triads, 

6) and to keep the code structure as close as possible to the scalar GCM, 


This vectorization of the Hydrodynamics has produced a remarkable overall 
increase in speed over the scalar code. The vector code streams through all 
the longitudinal and the vertical grid-points at the same time and, where 
appropriate, the latitudinal index as well. The vector Shapiro filtering 
process on the CYBER 205 is over 80 times faster then the scalar version on 
AMDAHL. Utilizing the U.K. Meteorological Office Fast Fourier- Transform code 
our Fourier filtering process gave us a factor of 25 over the scalar process 
on AMDAHL. We are also using a fast method to compute the expression, pressure 
to the kappa power. 

The vectorization of the hydrodynamic code is now complete and further 
improvements and enhancements of the model are currently being investigated. 
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VECTORIZEn SCM OBJECTIVE ANALYSIS SCHEME ON THE CYBER 205 
J, MoHni, W. E. Baker, 0, Edelmann, and J. Woollen 


The successive correction method (SCM) of objective analysis has been 
developed for assimilation studies with the FGGE data (Baker, 1983). In the 
SCM, successive modifications are made to the first guess fields provided by 
the fourth-order 6LAS GCM (Kalnay-Ri vas and Hoitsma, 1979), Eastward v<nd 
northward wind components, geopotential height, and relative humidity are 
analyzed on mandatory pressure levels. Surface pressure and temperature are 
reduced to sea level and analyzed there. The assimilation procedure involves 
the intermittent analysis of batches of data grouped in a h window about 
each synoptic time. 

In order to write efficient, vectorized’ code for the Cyber 205, the opera- 
tional code on the Amdahl V6 had to be substantially rewritten. The major 
internal data structures were redesigned. More of the working variables were 
saved in vectors. Time-consuming calculations, particularly the calculation 
of the gridpoint corrections, were rewritten using explicit vector syntax. The 
preprocessing and data checking were also modified to run as a separate program, 
thiis simplifying the I/O on the Cyber 205. Table 1 summarizes the CPU required 
for one analysis cycle. As may be seen in Table 1, executing the scalar (Amdahl) 
code on the Cyber resulted in only a modest improvement in performance. However, 
redesigning the computations to take advantage of the vector architecture 
resulted in a significant decrease in the required CPU time. 


Table 1. CPU required for one analysis cycle in sec. 


Amdahl Cyber Cyber 

(Scalar) (Vector) 

Surface analysis 

Height and moisture 
analysis 

Wind analysis 

72.4 44.5 10.2 

568.8 476.0 30.4 

388.9 206.3 28.4 

Total 

1030.1 726.8 69.0 
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June 22, 1982 


SECONDARY DISTURBANCES IN POLAR AIR STREAMS; 
COMMA CLOUDS . TROUGHS AMD POLAR LOWS 

R. J. Reed/University of Washington 


During the period October to April in the North Pacific comma- 
shaped cloud formations indicative of small synoptic-scale or large 
mesoscale disurbances occur frequently in polar air streams behind 
or poleward of major frontal bands. These systems are important to 
weather at sea and also at coastal localities when they come ashore. 
On surface weather maps they are generally marked by a trough that 
lies along the tail of the comma (depicted by a broken line and the 
symbol TROF on NMC analyses) and often by a low pressure center that 
lies beneath the comma head. On some occasions the trough acquires 
frontal characteristics and quite frequently the comma cloud excites 
and interacts with a wave on the polar front. 

In some cases, at least, the disturbances appear to resemble 
polar lows, small low pressure systems that move southward from the 
Greenland-Iceland area and bring showery weather to the British 
Isles, North Sea and adjacent continent. Also in some cases they 
may be associated with the bent-back occlusions and secondary cold 
fronts that at one time were commonly analyzed on synoptic charts. 

It is this host of phenomena, all characterized by development 
in polar air streams behind major frontal cyclones, that forms the 
topic of the present lecture. The main features of the systems will 
be reviewed and examples of several different types of developments 
will be shown. Finally, the origin of the systems will be discussed 
and pertinent theoretical and modeling studies mentioned. 
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July 13, 1982 


A SEMI-LAGRANGIAN AND SEMI-IMPLICIT NUMERICAL INTEGRATION SCHEME 
FOR THE PRIMITIVE METEOROLOGICAL EQUATIONS 

A. Robert/Canadian Meteorological Center 


A seml-lagrangian algorithm is associated with the semi-impl icit 
method in the integration of the shallow water equations on a rotating 
sphere. The resulting model is unconditionally stable and can be inte- 
grated with rather large time steps. Truncation errors remain reasonably 
small with time steps 25 times as large as those used with explicit 
integration schemes. 

An analysis of the proposed method is performed and it indicates 
that the scheme is stable. Also, the results of a few integrations are 
presented and from these we conclude that the model is not very sensUive 
to the size of the time step provided that it does not exceed a value of 
the order of two or three hours. 




July 20, 1982 


THE INITIAL SPECIFICATION OF THE ATMOSPHERIC STATS FOR 
NUMERICAL PREDICTION MODELS 

A. Kasahara/National Center for Atmospheric Research 


A deficiency in numerical weather prediction still exists in the gap 
between observational data analysis and the use of such an analysis as the 
initial condition for forecast models. In the traditional practice of 
analyzing meteorological data, one tacitly assumes that observational 
errors are often substantial and, therefore, observed data may be modified 
in the event that they don't agree with our models of the atmosphere for 
analysis and prediction. This practice was understandable in those days 
when meteorological observations were relatively primitive and scarce, but 
this viewpoint must change when dealing with more accurate and abundant 
observations than, say, ten years ago. 

The problem of initial specification of the atmospheric state for 
numerical prediction models will be reviewed. In recent years, notable 
progress has been made in our understanding of initialization to set up 
balanced initial conditions for numerical prediction models. Discussion 
will be centered on the relationship between two recent approaches, called 
the nonlinear normal mode procedure and the bounded derivative method in 
relation to the classical balancing based on quasi -geostrophic theory. 
Since both new approaches of initialization are more general than the 
classical procedures, the connection of the two approaches with the 
quasi-geostrophic formulation will shed light on closing the gap between 
data analysts and the preparation of initial conditions for prediction 
models. 



August 3, 1982 


ORIGIN OF MONSOON DEPRESSION 
R, Lindzen/Harvard University 


The development of monsoon depressions over the Bay of Bengal is a 
striking phenomenon restricted largely to the month of July. It is shown 
that such development cannot arise from either CISK or baroclinic insta- 
bility. Apart from the intrinsic weakness of baroclinic conversions, the 
momentum exchange associated with cumulus convection effectively inhibits 
both processes. This leaves barotropic instability as the likely mechanism. 

Studies of barotropic instability have thus far been restricted to 
normal mode analyses. Normal mode instabilities, however, are at variance 
with observed monsoon depressions. The most rapidly growing normal modes 
are associated with 200 mb winds, travel far too rapidly westward and do 
not penetrate into the lower troposphere sufficiently to organize convection; 
they, morecver, display no special response for July conditions. 

Instead of normal modes we examine the asymptotic response to pulse 
perturbations. Such asymptotics emphasize instabilities with zero group 
velocity. The calculation of such asymptotics is greatly facilitated by the 
existence of analytic expressions for stability relatione (phase speed as a 
function of wavenumber) in the complex plane. Fortunately, we find that a 
very sizable class of reasonable barotropic profiles yield stability problems 
homomorphic to the Charney problem in baroclinic instability. Analytic 
solutions to the latter problem have recently been developed with WKB methods. 
Such methods are used to show that instability to localized perturbations 
emphasizes the lower troposphere and is strongest in July, 


August 10i» 1982 


WHAT CAN SATELLITES CONTRIBUTE TO AN IfNDERSTANOlNa OF 
THE EL NINO PROBLEM? 

Klaus Wyrtki /University of Hawaii 


El Nino Is a sporadically occurring event In the ocean-atmosphere system 
of the Pacific Ocean and has catastrophic consequences for the fishery off 
Peru, causes torrential rainfalls In the otherwise dry coastal areas of 
South America and has repercussions for the climate over North America in 
successive seasons. The last major El Nino events were In 1957-58, 1965, 
1972-73, and 1976 and were apparently somewhat more frequent than in 
preceding decades. Research In recent years has been able to explain the 
oceans response by atmospheric forcing over the central equatorial Pacific 
and numerical models have succeeded to simulate the events for a prescribed 
wind forcing. Based on the dynamic theory of El Nino, which requires a 
build-up of warm water and higher sea level In the western Pacific caused 
by sustained strong southeast trade winds before the event, a prediction of 
the occurrence of El Nino becomes possible. Consequently, it will be 
possible to give advance warning of the next event, and this opportunity 
should be used to observe the entire sequence of atmospheric and oceanic 
events during the next El Nino and to obtain a set of observations never 
collected before. Satellites will play an Important role in monitoring the 
next El Nino event. 

Development of a strong El Nino event occurs In five phases. The build 
-up phase leading to El Nino Is characterized by a period of sustained 
abnormally strong southeast trade winds lasting more than a year. As a 
consequence, westward flowing equatorial currents are Intensified, leading 
to an accumulation of warm water in the western pacific and to higher sea 
level. At the same time upwelling along the equator and along the coast 
of Peru Is strong and the surface waters are cool. All these effects can 
be easily monitored by satellite winds, XBT sections from ships of opportu- 
nity and sea level stations. 

The triggering of El Nino through the collapse of the southeast trade 
wind field in the western and central equatorial Pacific usually occurs 
between December and March, during southern summer, as part of the seasonal 
relaxation of the wind. The extent and timing of this relaxation needs to 
be observed and studied, since it triggers the equatorial Kelvin wave, which 
travels east and causes downwelling along the coast of South America. While 
the deformation of the upper layer and of the discontinuity later In response 
to the Kelvinwave have been observed and modeled, the velocity distribution 
associated with It has never been measured and this would be a critical 
test for the validity of the model. The eastward advance of the Kelvin 
wave or group of Kelvin waves Is apparently also accompanied by a substantial 
decrease of the strength of the South Equatorial Current. 

About 6 to 8 weeks after Its triggering, the Kelvin wave impinges on 
South America, accumulating surface water along the coast and depressing the 
thermocllne. This documents Itself In a rapid and virtually simultaneous 
rise of sea level from Ecuador to Lima. The warmer coastal temperatures 
can be readily explained by the suppression of upwelling and the shoreward 
advance of warmer offshore water, but the quantitative aspects of the 
process need to be checked. Much more difficult to explain are the rising 
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surface temperatures farther offshore In the general area between the 
Galapagos Islands and about 15®S. The warming of these waters can be caused 
either by larger heat input at the sea surface (reduced cloud cover, less 
evaporation) or by a weaker flow in the South Equatorial Current causing 
less advection of collo water from the south and a longer residence time of 
the surface water in the area, which would lead to in situ warming. Another 
aspect of El Nioo during this phase is the overflow of warm low salinity 
water from north of the Galapagos Front towards the south between the coast 
and the Galapagos Islands. This overflow was observed in 1975 and is 
probably a feature occurring annually on a small scale and during El Nino on 
a very large scale. 

The third phase of El Nino is characterized by a poleward spreading of 
the thermocline disturbance, by a westward spreading of the sea surface 
temperature anomaly along the equator, by a continuing draining of warm 
water from the western Pacific and lowering of the sea level, and by 
intensified eastward flow especially in the North Equatorial Countercurrent. 
Surface temperature anomalies appear first off the coast of Peru and spread 
westward reaching the Line Islands about five months later. It is not very 
likely that this spreading is an effect of advection, but more likely that 
the suppression of equatorial upwelling progresses westward as sea level 
in the equatorial trough rises. The time variatiion of equatorial upwelling 
and in particular the development of east-west slope of thermocline topography 
and sea level along the equator needs study- 

The continuous draining of water frort*; the western Pacific has been docu- 
mented well by both sea level measurements and by the response of the thermal 
structure observed from XBT sections and the 137* hydrographic section in 
the western Pacific. The response of the equatorial current system in the 
central Pacific has been monitored by means of sea level gauges. It consists 
of a weakening of the westward flowing South Equatorial Current and an 
intensification of the eastward flowing North Equatorial Countercurrent and 
probably also of the South Equatorial Countercurrent. Little is known 
about the response of the Equatorial Undercurrent, which will require direct 
observation. During this phase of El Nino, which corresponds to southern 
winter, the sea surface temperatures off South America drop again as upwelling 
resumes under the influence of seasonally stronger coastal winds. A weak or 
moderately strong El Nino might terminate during this phase. 

During the fourth phase of a strong El Nino event the eastward flowing 
North Equatorial Countercurrent advects much warm water into the Eastern 
Tropical Pacific, leading to considerable warming there. This strong flow 
together with a reflected wave from the western boundary may be the reason 
for the second peak in sea level, thermocline depression and sea surface 
temperature observed during strong El Nino events. 

The fifth and last phase of El Nino is the unexpectedly sudden termina- 
tion of the event and a rapid return to normal conditions. Sea level drops 
along the coast of South America, rises even more rapidly in the western 
Pacific, equatorial upwelling resumes, surface temperatures become cooler 
than normal. All happens within a rather short interval of about two 
months. The processes involved in this rapid return to normal have not 
really been investigated. 
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Satellites are expected to contribute significantly to the monitoring of 
the next El Ml no event. Most Important Is the determination of the wind 
field over the ocean, which can be accomplished by either determination of 
cloud motion vectors or by scatterometer measurements of wind stress. 

Improved methods to observe sea surface temperature will allow a better 
determination of this Important field, especially when appropriate corrections 
for water vapor can be applied. Water movements can be observed by the use of 
satellite-tracked drifting buoys, which might also relay other environmental 
data. Satellite observations "of atmospheric circulation during and before 
an El Nino event are of course needed to Interpret the oceanographic data 
and to study the Interaction between the two systems. 

On the other hand, measurements of subsurface thermal structure are not 
yet possible from satallltes, and direct methods will have to be used. This 
applies as well for the collection of most chemical and biological data. 

The greatest advantage of the satellite Is, however, the fact that It 
delivers a contlnous spaclal coverage of global extent. This fact alone 
will revolutionize oceanographic research and thinking In the years to come. 
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